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I . INT R ODUCTI ON 



The object of this thesis has been twofold. 

The first object was to develop FCSTEAN versions of 
several existing APL programs which were designed to analyze 
univariate data, in particular to test fez exponentiality 
and normality of the data and, by sectioning or jackrifing 
obtain estimates of variances of sample moments. 



The second object was to use these programs in a 
simulation study of first-come first-served gueues in which 
the service times and the inter-arrival times are 
exponentially distributed but dependent. The dependence is 
introduced by using the mixed moving average autoregressive 
structure ( EARMA(p,g) ) for exponential sequences 
introduced by Gaver and Lewis (1977), Lawrance and Lewis 
(1977), Jacobs and Lewis (1977) and Lawrence and Lewis 
(1978). four schemes are introduced which give 
autocorr elated and cross-ccrrela tad service and arrival 
times in various degrees. 



Since nc analytical properties or these queues can be 
derived, their properties are studied by simulation. 
Because the EAEMA gueues are neither regenerative nor 
ilarkovian there are several problems in the simulation study 
of knowing when one can assume that a steady state has been 
reached in the simulation and of obtaining confidence 
intervals for estimated parameters. This problem is handled 
by generating a number of independent sample paths and 
comparing statistics obtained from the replicatiers at 
several points along the sample paths. Only the waiting 
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not the state 



time process (W ) is studied in this thesis, 

N 

space process N (t) . 

Unfortunately the simulation study cf the EAEMJ-type 
queue was hampered by rhe size and speed cf the 
Ian System 3c0/67 system at the Naval Postgraduate School. 
Consequently the results were not as extensive as had been 
hoped fcr. Fairly detailed results are given for the case 
where the service time process is autoregressive ; one case 
where the traffic intensity is 0.99 and the correlation 
parameter 0.98, was intensively investigated. This gives 
seme idea of the length of the transient in the queue and of 
the inflation of the mean waiting time because cf the 
correlation in the service process. Some investigation of 
the case where the service times and inter-arrival times are 
cross-correlated are given. 

In section II we describe the basic histogram and 
statistics computation package 'HISIGS/F' and apply it to 
two sets of telephone error data. In section III the 
plotting method and the plotting Subroutines 'NGFMEL', 
•HXPLT* are described and we test them with various 
generated data. In section IV we describe the assessment of 
variability Subroutines 'S2CTN', 'JACK' and we apply them to 
telephone data set 1. In section V the Subreutine 'LIST' is 
described and we apply it to two sets of telephone error 
data. In section VI we describe and analyze rhe EAEMA(p,g) 
model as well the program which we developed in order to 
simulare it. Computer program listings are provided after 
section VI. 
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II 



SUBROUTINE HISTGS/HISIES 



A. DESCRIPTION 



Here 


the 


Subtoutine 'HISTGS/HISTFS' 


is 


press 


nted. 


'HISTGS/ 


HISIFS' 


is used for estimating 


the 


probab 


ility 


density 


function 


from a given set of data 


and 


at the 


same 


time CO 


mputing 


some basic statistics. 


Basically 


this 



Subroutine is the library FORTRAN-Subroutine HISTF/G, which 
was created at N.P.S. by B. H. RobinsoE. It has been 
modified by adding the new entry point 'LIMITS'. 

The new entry poinx may be used for wild data or for 
data having a mixed distributicn. Therefore it is a useful 
tool for that user involved with this kind of data. Simply, 
the user has to define those values in the range of data 
(sectioning the data into several disjoint sections) that 
he/she believes are useful cut points. For each section of 
data the user may have the probability density function as 
well as the basic statistics ( as HISTF/G does ) . Also by 
sectioning the data over its range, the user may have a 
histogram of the random variable when it is conditioned to 
be between given limits. 



The number of 
accept is as many 



sections that the entry point 'LIMITS' 
as 5 C . 



can 



A complete 
is given in the 
the subroutine 



description of 
subroutine. On 
is given by 



how ' HISTGS/HISTFS ' operates 
the other hand, a summary of 
typing on the terminal the 
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command CESCEIBE HISTGS under CMS (Camtridce Monitor Sjstsm, 
Hef, 4). Typing 'DESCRIBE HISTGS' the following response 
will be given on the terminal: 

SUBfiCOTIKE HISTGS/HISTES 

' HIS IGS/HISTFS ' gives you a histogram of a set of data 
along with the estimate of a set of basic statistics. also 
it gives you the following options: 

1. Displaying a smoothed empirical probability density 
function over the histogram, 

2. Sectioning the data into several disjoint sections 
(no more than 50) and then having a hisrcgram and the 
estimate of basic sratistics for each section, 

3. Scaling the histogram, 

4. Displaying just the histogram (with or without the 
density function) and no statistics. 



The estimated basic statistics are: Mean, M 

Irimean, Midrange, Geometric Mean, Harmonic Mean, Var 
Standard Deviaricn, Coefficient cf Variation, 
Deviation, Range, Midspread, Third and Forth Central M 
Skewness, Kurtosis, Beta 1 , Beta2, Maximum, Minimu 
Quantiles . 



edian , 
iance. 
Mean 
cment , 
n , and 



CALLING 3 


SQUENCES 




CALL 


HISTGS 


( X, N, 


NBAR 


CALL 


HISIES 


( X, N, 


NBAR 


CALI 


LIMITS 


( X, N, 


XLIM 


CALL 


NCSTAS 






CALL 


STATS 






CALL 


FIIS ( 


SCALE ) 




CALL 


NCFIXS 






CALL 


PECM AX 


( PSC ) 




CALL 


NCERMX 







19 



ARGUHENTS 



X 


Array of 


data values 


N 


Number of 


data values 


NEAR 


Number of 


bars in the histogram 


XLIM 


Array of 


cut points 


LIMS- 


Number of 


cut points 


SCALE 


Vector of 


two values to scale the 



PSC A real variable (between 0. and 1.) defining 
the inaxiiaum value of the probability axis. 

If the entry point NOSTAS is called before callinc any 
of the entry points HISTFS, HISTGS, or LIMITS then no 
sxatistics are printed, otherwise xhe statistics are printed 
by default. 

By calling the entry point FIXS the histcgram scale nay 
be fixed and remains set unless it is reset ty another call 
to FIXS or allowed to float (the default) by calling NCFIXS. 
The reason for fixing the scales is so that when ccmparing 
mere than one batch of data, we can get ccaparably scaled 
histograms. 

By calling PBOMAX, the maximum value of the probahility 
axis (Y) can be set to PSC, which will remain in effect 
until NOPRMX is called (the default) . 

More information is given in tne sucroutine. 

E. USING ' HISTGS/HISTFS* WITH TELEPHONE DATA 1 AND 
TELEPHONE DATA 2 



Here 'HISTGS/HISTFS* has been used to plot the set of 
Telephone Data 1 and Telephone Data 2. The algorithmic 
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the 



procedure for both sets of data is the same. That is, 
entry point *HISTFS* is called first and then the entry 
point 'LIMITS*. Four cut points have been used (1, 2, 141, 

86000) sectioning the data into three disjoint sections as 
fellows; 

Section 1 : Consists of data points X , such that 

< i> 

1 < X <2. 

<i> 

Section 2 ; Includes the data pcinrs X , such that 

<i) 

2 < X < 141. 

<i> 

Section 3 : Includes the data points X , such that 

<i) 

142 < X < 86000. 

Ci) 

No print-out for section 1 is given because cf the 
constant value of the data points. 

Figures 1 through 6 give the histogram and the basic 
statistics fer both sets of data ( For rhe entire set and 
fer sections 2 and 3 ) . 

Me observe that both of data sets are sc unruly that the 
use of the entry pcint 'LIMITS' cf 'HISTGS' is suggested. 

Observing Figures 1 and 4 we may conclude that both data 
sets appear to he the same, having seme geometric or 
exponential distribution. But as we will see later this is 
not true. Thus the histogram, in that case, which is 
obtained by using just the entry point 'HISTFS' is not 
helpful and leads us tc a wrong conclusion. 

'LIMITS' is applied now to give us a useful answer : 
Observing figures 2 and 5, obtained by using the entry 
point 'LIMITS', we may conclude that the data sets dc not 
have the same distribution ( at least in that interval ). 
In addition data set 2 is far from having a geometric cr an 
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seen 



exponential distribution, as it can be seen by the histogram 
of the entire set of data. This is a conclusion which 
cannot be obtained by without use of the entry point 
•LIMITS* . 

Observing also Figures 3 and 6 , where the histogram of 
the 3rd section from each data set is presented, we may 
conclude that this section of both data sets may have an 
exponential distribution, but not the same one. 

The mosx important facet of the data which appears is 
the modes at anout 124. The mode is much more apparent in 
Telephone Dana 2. In fact closer investigation reveals that 
this mode was due to leakage of dial pulses; this leakage 
caused bit errors at multiples of 124 bits. 

The above is a quick and informal analysis of the 
Telephone Data Sets obtained by using the entry point 
•LIMITS' of • HISTGS/HISTFS* , and the conclusion is that both 
of data sets are so unruly that any formal oharacteriration 
cf their distribution is not so easily obtained. In 
addition it is apparent that a mixed distribution governs 
both sets of data. 



Having quickly analyzed both of the telephone data sets, 
the use of entry point 'LIMITS' cf the Subroutine 
' aiSTGS/aiSTFS' is obvious. 
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ffiEQOtNCIES SAUPLE SIZE = 672 
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FBECUtNCIES SABPLE SIZE » 406 
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F 8 ECUEMCIES SAflPLE SIZE = 138 
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PBECUfHCIIS SiHPLE SIZE - 736 
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FBEgUENCIES SAMPLE SIZE = 315 
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FBEOUtNCIFS SAMPLE SIZE * 243 
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Ill 



PLOTIING SUBROUIINES 



k. INTBOCUCIION 



Graphical methods of assessing the validity of a 
probability model, and of estimating some basic parameters 
(especially location and scale parameters) have been ^*idely 
applied over many years. A variety of different 
prescriptions have been advanced for the plotting positions, 
that is, fcr the set of values at which the ordered 
cbservaticns in a sample should he plotted. 



Plotting 
visual methc 
distribution 
about the d 
probability p 



the ordered data 


is an in 


formal and 


quick 


d for getting an 


impression 


of 


how 


the 


of data looks. 


Depending 


cn what is 


known 


istributicn of the 


data xh 


ere 


are 


many 



lotting methods. 



The 

•EXPLT* 



Plotting Subroutines presented 
use the method cf plotting 



here, 'NORflPL' and 
the order statistics 



( X ) versus their expected value ( EfX ] ) called 
<i) Ci) 

•SCORES'. This method does not require knowledge cf the 

Cumulative Distribution Function F (x) (except for its 

assumed continuity) . Therefore it may he very attractive 

fcr those who knew nothing about the data to he analyzed. In 

addition the indication of the graph may lead the analyst to 

further useful results. 
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Using the Plotting Subroutines 'NORMPL' and 'EXPLI' one 
may analyze the resulting graphs obtained tc get an idea of 
the distribution of the data. 



E. DESCBIPTICN OF THE PLOTTIHG METHOD 



The idea behind this method, as above stated, is to plot 

X versus £[ X ]. 

Ci) (i)' 

It is used where data X , X , ..., X arise as 

12 n 

independent observations or a continuous Bacdom Variable X 

with a distribution function which is believed to have some 
particular form F (z) . Then the ordered data 

X , X , X are plotted against Efx ]. 

<1) <2> tn) <i) 

If the model for F(x) used to get E[ X 1 is true then 

( i ) 

the resulting plot is a linear (regression) relationship and 

should plot, within limits of sample variations, as a 

straight line. Barnett [2] gives the value of E( X ] as a 

^ i J 

linear relation of the plotting position That is: 

E[X ] = m + sx , 

here m, and s, are the location and scale parameters. The 
Mean and Variance of X are not necessary. As we see the 

above relation uniquely defines the plotting position as : 
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Ef X -ml 

(i) 

X = 

<i) s 



X - m 

= E[ ] , 

s 

Therefore, if the model is correct the linearity will 
frobahly be evident. In addition if the resulting plot not 
being linear leads us to reject the model, we may cet an 
indication of the type cf departure present and thus an 
indication cf what action to take or what alternative model 
to use. It should be noted that the linearity of the graph 
suggests informally that the model is true, but not to 
accept the model. Additional formal tests must be done in 
order to accept the model. 

C. SUEHCOTINE NCRMEl 



1 . Cescriptio n 



Subroutine 'MCRMPL 
independent observations of 
having a distribution funct 
Normal form. 



, is used to test 
a continuous Random 
on which is believed 



a sample of 
Variable X 
to have the 







• NORIfEL' 


sorts the 


data into 


increasing 


cr 


cer to 


create 


a 


set cf 


Order Statis 


tics. Then 


the Order S 


zd.t 


istics 


are 




plotted 


versus 


either 


Normal 




Scores 


(if n 


< 


50, n 


= number of 


data-points) 


or inverse 


0 


i the 


Standard 


Normal 


C umulati ve 


Distributi on 


Function 


(F 


1 

(y) ) / 


(if n 


> 


50) . 
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Evaluation of Normal Scores is based on computing 

t h 

the Expected Value of the i Order Statistic ( E[X^,^] ) . 



This evaluation follows: 

Let 

t h 

g(x^.^) be the density function of the i Order 



Statistic of a sample 1 , X 1 from a population 

1 2 n 



having the Standard Normal Cistribution F. Then g(x^_^) is 



given by : 



n< i-1 n-i 

ci)' (n-i)i X <i) X (i) X (i) 

n» i-1 n-i 

P Q Z , 

(i-1)| (n-i)! 

where : 

2 

-X /2 

1 (i> 

Z = f (X ) a , 

X <i) 1/2 

<2 R) 

u; 2 

1 -u /2 

P = F (X ) = e du , 

X Ci> 1/2 

^00(2 n ) 

Q = 1-P . 

Now 

E£ X^i^ 3 = j X g (X) dx 

-Oo 

n' i-1 n-i 

(1) = 2 X Z P Q dx . 

(i-1); (n-i): J 

^-co 
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The value of expression (1) has been evaluated and 
is given in [19]. This table gives the values of Normal 
Scores frcm 1 to n/2 (for n<50) . The values from (n/2) + 1 to 
n are evaluated using the symmetric expression : 



Also 



E[X ] = 
Ci>-' 



■ ■ 1 J 

<■ n-i+ 1 J 



E[X ] = 0, n = 2k+1 , 

(Jc+1) 



•NOEiUFL* uses the above Normal scores if n < 50. For 
n > 50 ‘NOfiaEL’ computes its cwn Normal Scores using the 
program function 'INyMSM', which gives the value x of the 
inverse Standard Normal distribution (x = F (y ) , 
0 < y < 1) . Justification of this is the fact that if the 
data are really normally distributed then 



E[ X ] F( i/ (n+ 1) ). Therefore 'NORMPL* can acconcdate 

t i J 

any sample sise. 



In addition to the plot given by 'NOEMPL', the 
Hilk-Shapiro test for normality can be evaluated by the same 
Subroutine. A derailed description of that test is gi^ien in 
[22]. However a summary of that test is presented here in 
order to help us to describe how the program works. The 
Silk-Shapiro test is based on the statistic: 

2 

b 




S 



where : 

k 

b = 7" a (X - X ) , k = c/2, n even 

i = 1 n-i+ 1 (n-i + 1 ) (i) 
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b = a (X -X ) <• a (X -X ) + . . . 

n <a> <1) n-1 <n-U <2> 



+ a (X -X ) , 

k + 2 <K + 2> <k> 



k= (n- 1) /2 , n odd 



2 n 2 

s = H (X - I) 
i= 1 i 



Ih€ coefficients (a ) are defined by; 

i 



n 

a = ZIau /C, i=1,2,...,n 

i j=1 j ij 

where ; 



m = E[X ], (That is the 'NORMAL SCORES') 

j ^ 3 ^ 

u = Ccv[X ,X ], i,j = 1,2, ...,n 

i j f i ) ( j) 






C = 



I _i _-j 

*a V V m 



where ; 



ID “ #1^ ) 

12 n 



V = (u ) , That is the nXn Covariance Matrix, 
ir j 



There are some approximations associated with the 
evaluation of these coefficients. However fer n < 5C the 



(a) are tabulated in the above-mentioned paper and the 



program uses them by FCHTRAN DATA initialization stateaents. 
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fcr n>50 the (a ) are computed by the program using 

i 

the following approximation method: 



a = 2m /C, i = , 

i i 

where C is given by the following Least-Square equation as a 
function of n: 

2 

C = -2.722 + 4.083n 



Also, for n=1 



2 + - 5 ) 

a 

1 1/2 

2 P(n/2 + 1) 



Note that a = a 

i n-i+1 



Some properties for the W Statistic can be given 
here, which have been taken from the above mentioned paper; 

(a) . ri is scale and origin invariant. 

(b) . has a distribution which depends only on the 

sample size n, for samples from a Normal Distribution. 

2 

(c) . w is statistically independent of S and of 1 for 



samples 


from 


a Normal Distribution. 




(d) . 


r 

E[M 


2r 2r 

] = E[ b ]/E[ S ] for any 


r . 


(e) • 


I he 


maximum value 


of w is 1 . 




(f) . 


The 


minimum value 


2 

of W is na 

1 


/<n-1) 


(g) • 


The 


half and first 


moments of 


W are 






1/2 

E[ii ] = 


2 

R Rn/2 - 


.5) 






C r(n/2) 


V2 

2 






T1 r TJ T — 


2 2 
R (R +1) 








EL W J — 


2 


• 



2 

C (n-1) 



given 






by ; 
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where; 

2 ' -1 

B = IB 7 m 

(h) . For n=3 the density of H is given as: 

3 - 1/2 - 1/2 

f (w) = --- (1-w) M , .75 < w < 1. 

W *' 

2. Program Structure 

•HOEIUPL' is a FOETH AN-callable Subroutine with each 
call returning (optionally) either a plot fcr a given set of 
data or a value for the W Statistic or both. 

The program is divided into three parts: 

The first part is the control prcgraa of 'NOFMEL'. 
In that part the 'Normal Scores* for n = 2(1)25 and 
n = 26(2)50 are stored along with the Subroutine by DATA 
statements; on the other hand, the Normal Scores for r > 50 
are computed by calling the pregram's function 'INVNEM'. In 
this part the user specifies 'ust a plot cr just the value 
of W or both. The Library-Subroutine 'PXSCBI' is used by 
the program. 

The second part of the program censisrs of the 
Subroutine 'PLOT' which accepts any data set to be tested 
fcr Normality. 'PLOT' itself scales the data according to 
its range and plots the scaled data along 110 equal spaced 
positions of the X-axis. No plot is given if the data has a 
constant value. 

The third part of the program is the 'WILK' 

Subroutine. In this Subroutine the coefficients (a ) 

i 

(for n < 50) of the W Statistic are stored using FCSTEAN 
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DATA initialization statements. For n>50 the coefficients 
(a ) are computed by the Subroutine itself. Ihis Subroutine 

also calls the 'INVNRM* Function to get the ’Normal Sccres* 

to be used fcr the ccmputaticn cf the (a ) *s vihen n > 50. 

i 

A ccmpiete description of how 'NOEllEL* operates is 
given in the Subroutine. However a summary is obtained by 
typing on the terminal DESCRIBE NOHUPL. When the user types 
the command DESCRIBE NCRMPL under the CMS environment the 
following response is printed on the terminal; 

SU3ECUTINE NCRMPL 



’NOEMPL' takes a set of data, sorts it inxo 
ascending order and uses the created Order Statistics for: 

1. Plotxing X versus either 'Normal Scores' (if 

(i) 

n < 50, n = Sample size) or Inverse cf the unit Ncrmal 
Cumulative Distribution (that is; 



-1 

X versus Normal scores or F ( i/(n + 1) ) ,) 

(i) 



to see if there is 


a linear fix. 








2. Computing 


the value of the 


w 


t est 


Statistic 


(Wilk-Shapirc test 


for normality) . 









It is called by; 

CALI NCEMPI ( X, SCCRES, N, K ) 



X 

SCORES 

N 

K 



ARGUMENTS 



Is the 


array co 


Is a work array 


Is the 


number o 


U ser ' s 


option 


For: 


K = 1 



ntaining the daxa 
of dimension N 
f data values 

A plot only , 
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K = 2 The H-value is given , 

K = 3 Both of the above are given. 

Mere information is given in the subroutine. 



3 • Interp ret ing The Output 

If data are really Normally distributed then a 
straight line is expected from the output. This is an 
indication that gives the user a first feeling of the 
distribution of data. But if the plot is not linear then 
the user has to reject the normality assumption. In 
addition tne shape of the plot may lead the user to make an 
alternative decision of what model he may use. For example: 

(a) . Shape of the form of Figure 7a suggests that 
data are skewed to one side (specifically positive values) . 
In that case an Exponential model could he an appropriate 
one. The histogram of the data is given in Figure 7b. 

(b) . Shape of the form of Figure 3a can be 
interpreted as data not as dispersed as the Normal 
Distribution. This set of data may have a Unifcrii or a 
Triangular D istriburion, and could, for instance, arise if 
measuring tolerances ox components hand-picked to lie within 
tolerances. The associated histogram is given in Figure 3b. 



(c) . Shape 
mere dispersed data 
probably coming 
Distribution. The 



of the form of Figure 9a indicates much 
than Normal, but the data is symmetric, 
from a Double Exponential or Cauchy 
associated histogram is given in Figure 



9b . 



net 



It should be noted here that a straight 
always arise from Normal data. Figure 10a 



line does 
for example 
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from a T 
ue obse 
results 
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of appl 
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indication that data are Normal 
£6 data have been generated by 
c Triangular Distribution in the ra 
er hand figure 11a plots a data set 
rianqular (0 to 2) distribution. E 
rve that the graph is not a stra 
demonstrate the informality of 
Method and in these cases the use 
the data by other formal tests. Ne 
thod is primarily an initial scree 
ication is important. Figure 
ed histogram. 



(Straight lice) . 

computer from a 
nge (0 to 2 ) . On 
cf size 50 again 
ut in this case 
igtt line. These 
that Probability 
r has to ccntinue 
vertheless, since 
Ding device, ease 
11b gives the 
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MUafiEB Of OBDERCD FAIBS - 200 



09 



Figure 7a - GENERATED DATA FROM AN EXPONENTIAL 

DISTRIBUTION 
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NUn3£B OF OBDEBED PAIRS « 200 



o 



o 
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Figure 8a - 



GENERATED DATA FROE A UNIFORM (0,1) 
DISTRIBUTION 
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NUfiBEa or OBDERCD FAIRS - 200 
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o 




Figure 9a - GENERATED BATA FROM A CAUCEY DISTRIBUTION 



4U 



-SCALE : = 0.859B-01 UNITS 



fBECOlBCIES SAMPLE SIZE * 200 



o 

<N 




03 



OOOOOOO 
I I 

Ck3MUj:>Juia}Cd 

O O ("n vO irt 
00004 
mp'O" r^r^C'® 



^ roC'^N CM 
lift 




O'-'O 

X.OX 

M Wl-4 



Z UCx)CdU<U 

O 

H He-t-e-6- 

3 ZZZZZ 

J3 C -< «5 -S -= -S 3C 

•M 3333333 

z x: cnoooox: 



fN 1 
(B t 
I 



tA zo*novOO>^ 
M -4*-o4vnr-c'-< 

a x: . « • . . s: 



CKzrzmczrz; zrzrxizrxrcxizrzrzrr; cz;z;z:z:x:z;z] i;z]z;x:x;x:z;z:3:x:z;z:z;x;z;z;z;z;z;£x:x;® 






in <n4s^oo<n3- 
H 300000 
z 

uucotuiUDa 

C •3’3‘Nfi®sO<7' 
O nOC'—'^nC® 
«z — O'l^uOrsiO 

•-r^^r^vrivn 

oo'mmvor* 

c :?• oioir^^rsj 



H 04^«-(T><N*- 

Z I I I 

aa 

o (n<n 

in>~* 

aa ct3in 

M zo*"<N 

s 

O COZg-'^ 

X x: z:cn:K:3a3 



I 



«r» 

rsi 



o 

04 





* * 




r-o*-0*-0 

OOOOOO 

M U bJ bJ U uJ 
\45 ^u^»/~o»r> 
<Nr*o siSOu^ 
OjO'OVrt ^sC 
'N<7'o<7'r^r' 

»“000»04i— 



(^U*»»-<Ni/3rs| 



3 

W m> < 
Uxeo bj 
a zu;>a z 

C «S 3 CO 0« 
bJ i-t o<zcn 
a zcoKzo 

co :»<nox«x 



*“<Nr4rMO 
OOOOO 
tilt 
UUiMIBbJ 
if ojmcT' 
mrsi004 

3 ^ 04 .T •-O' 

u r'r'0>«r>*“ 

z (MvyrMr^ix 

M 

3 iDonf" if tO 
Z I I I I I 
M 



zzu 
«« z ^ ^z 

Z <<bJbJ< 

H z*-*zx:z 

Z <Ui-iUO 
M 

ij XZ4r>ZC 



Figure 9b - 



HISTOGRAM OF THE DATA Cl FIGURE 9a 



45 



NQHBEB OE OBDEBED FAIBS « 200 
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Figure 10b - HISTOGRAM OF THE DATA Cf FIGQRE 10a 
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Figure 11a - GENERATED DATA FROM A TRIANGULAR (C,1) 
DISTRIBUTION. SAMPLE SIZE 50 
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WILK-SHAPIHO TEST: TEST VALUE « 9.6880-01 SEE BIOttETRlKA (1965), V. 52, P. 591, FOB CBITICAL VALUES 



FBECUtHCIES SAMPLE SIZE 
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HISTOGRAM OF THE 
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4 



• * NOEMPL ' with Generate d Data fro m Various 
Distributions . 

The Subroutine ’NOfiaiL* has been used here to test 
for Noraality some data sets which were generated ty the 
ocmputer froa various distributions. The plots of these 
data sets illustrate the advantages and disadvantages of the 
plotting Subroutine 'NCBMPL'. 

Observing Figures 7a-b, 8a-b, 9a-b and 12a-b we may 
note the following: 

(a) . Figure 12a shows a plot of a sample of 5C data 
points from a Normal Distribution ( N (1000,1) ). As it was 
expected a very straight line is fitted. In addition to 
that plot, the Subroutine 'NORMPL* uses the option to 
evaluate the ' WILK-SHAEIRO' test and the W-value is cciiputed 
and printed on the same figure. Comparing this value with 
the percentage points of the W-tes~ which are given in [22], 
for sample size n=50 we see that we may accep* the normality 
assumption with a significanse level = .02 . Figure 12b 
gives the associated histogram. 

(b) . Figure 7a plots a sample of 200 data points 
of the Unit Exponential distribution. Here the nonlinear 
graph is obvious ( as is expected ) and suggests to 
reject the normality assumption of the data without doing 
any further formal test. 

(c) . Figure 8a gives a plot of a sample cf size 
200 from Uniform ( C, 1 ) generated data. The shape of 
that plot indicates data not as dispersed as the normal, 
therefore the user is told to reject normality, getting at 
the same time a suggestion that data may be Uniformly 
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distributed ( because of that particular shap 
graph ) . 



(d) . Figure 9a suggests str 
Normality assumption. Indeed this set 
generated from a Cauchy Distributio 
indicates a symmetric distribution. 
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NUBDER or ORDERED PAIRS 
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Figure 12a - 



GENERATED DATA FROM A 
DISTRI30TI0N 



NORMAL 



(100C, 1) 
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RILK-SHAPIRO TEST: TEST VALUE * 9.902l>-01 SEE BIOflETRIKA (1965), V. 52, P. 591, EOR CRITICAL VALUES 
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Figure 12b - HISrOGRAM OF THE DATA OE FIGURE 12a 
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5 



• ISORJPil With 



Cost Overruns Eata 



was used with the cost overruns data 
to see if there is a linear fit. Leo 



'NORMPL' 

for the year 1950 
at Figure 13 we 
since the plot is 
points deviate fr 
be conside 



see that data may be 
very close to linear, 
cm the straight line. 



[ 6 ], 

king 

Normally distributed 
Of course seme data 
but these points may 



red as outliers. 
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NUn^KB Qt ORDEBED FAIRS 



Figure 13 - 



COST OVERRUNS DATA 
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-SCALE : ^ 0.637E-01 UNITS 



D. 



SUBaCUTINE EXPIT 



1 . Descri ptio n 

Subroutine *EXPLT' is indended tc be usee for 

testing a saaple of independent observations of a continuous 

random variable X with a distribution function which is 

believed to be Exponential, by plotting each of the Order 

Statistics X versus its Expected Value, and by estinating 
U) 

the parameters Jj^and of the observatiens ( and 

introduced later) . 

This Subroutine sorts the data into increasing 
order, obtaining the Order Statistics. The Order Statistics 
are then plotted versus Exponential Scores. These scores are 
evaluated by the Subroutine. The Exponential Scores are the 
expected value of the Order Statistics ( E[X ] )• They 

(i) 

are more easily derived and computed than the Normal Scores. 
That is: 

Let X be a random variable Exponentially distributed 
with parameter > = 1 
Then it is known; 

f (X) = X > 0 , 

X 

F (X) = 1-e " , X > 0. 

X 
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th 

The Density c£ the i Order Statistic is given by: 

n-1 i-1 n-i 

f (X ) = n( ) F (X ) (1-F (X )) f (X ) 

X ' Ci)' 'i-1 X ti> X <i> X <i) 

<i> 

St 

Ibus for the 1 Order Statistic we have: 

- jx -;pc 

( 1 ) a- 1 it) 

f.. J = n(1- (1-e ) ) ys 



( 1 ) 



i 1 ) 



-n>x 



= nxe 



( 1 ) 



( 1 ) 



Therefore X 

f 

parameter "> = n . 
And for ':x= 1 have: 



is exponentially distributed with 



E[ X ] 
( 1 ) 



1 _ 1 
n \ n 



It can be proved (see Feller [7]) that the n random 

variables X , X -X , X -X , . . . , 1 -X from 

(1) (2) (1) (3) (2) (n) (n-1) 

an Exponential random variable X with densitj function: 

f (X) = >e , X > 0 , 

X 

are independent and the density of the random variable 

1-1 - X is given by; 

(k) (k-1) 



- :\(n-k+1) y 

f^(y) = (n-k + 1) :^e 



y > 0 



That is, Y is an Exponential random variable with 
parameter (r-k+1) TV . 
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Iberefore 2[ Y ] is given by: 



1 

E[ X -X ] 

(k) (k-1) :^(n-k + 1) 



1 

n-k+1 ' 



for =1 



Ibis fact is used to derive E[X ], i = 2,3,...,n . 

(i) 

I hat is: 



E[X ) = E[ (X -X 



( 1 ) 



( 1 ) ( 1 - 1 ) 



) + (X 



(i-1) ^(i-2) 



) + . . . + 



(X -X ) + X ] 

( 2 ) ( 1 ) ( 1 ) 



Ef X ] + E[ X -X ] + ...+ 

( 1 ) ( 2 ) ( 1 ) 

E[X -X ] + E[ X -X ] 

(i-1) (i-2) (i) (i-1) 



111 1 
n n-2+1 n-3+1 n-i+1 



1 

k=1 n-k+1 



i - 2 , . . . n 



'EXPLT' derives the values of E[ X ] using the 

<i) 

anove formula and then it uses its Subroutine 'EPLCI' to 



plot the Order Statistics versus their expected value. 



In addition to the plot, which is an informal and 
quick test for exponen tiality the estimates of the 
statistics and ^ 2 . Provide another informal test since the 
values of and have a constant value for any 
Exponential distribution: 
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■ ' i ’ 

• i i! 



J.;it E 

sfiy. E 



I 




Ey definitioE ^j_and are given by 



Ih us 



<7 



- 3 . 



EC 

cr^ 



1 3 2 2 3 

-- E[X - 3X + 3X y. - ] 



(1) 



1 3 2 2 3 

-- (E[X ] - 3yE[X ] 3 y E[X] -y ) 



Now 



E[X ] = 



3 

X f (X) dx 
<7 X 



3 -:\x 

= 3o ^ 



r(4) 



^3 



dx 

4-1 4 -;\x 

X e 

r(4) 



dx 



Ihe above integral is that of the Gamma density fun( 
with parameter OL = H , and since its value equals 1 
have : 



3 31 6 

E[X ] = — = 



tion 

we 



59 



Ccaiputing E[X ] we have; 



E[X ] = 



« 2 


-;\x 




e dx 


^ O 






3-1 


r (3) 


r.f— 


■•'i" 


i r(3) 



^8 dx 



Also 



0^^ 

Z[X] = 4- = 



Cv V 

Substituting in (1) we have: 



0 (x-'r*) 

h = EC— 1 



_ 1 6 12 11 






3 3 

= Ov { 2/^ ) = 2 . 

Iherefore the value of "ij ^ equals 2 for any Exponential 
distribuxicE . 

Also is computing as follows: 

i - 3 

4 3 2 2 3 4 4 

= E[ X - 4y»X +6yX - 4yX+y]/cr - 3 



( 2 ) 



4 3 2 2 3 4 4 

{E[X ] - 4yE[X 3 + 6y E[X ] - 4y E[X] + y } /<:t-3. 
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£at 



4 (CO 4 -;\x 

j = \ X dx 



4 



oo 5-1 5 -:^X 

X ye dx 



r(5) 



= 21^ / y 



Ihus 



24 



^1= "7 ^ — 






12 11 1 
+ 6 — 4 + } 



y 









- 3 



4 4 

= 1 (9/^)-3 = 6 



Therefore the value of equals 6 for any 

Exponential c istr ibuticn. 

Ibis fact is used by the Subroutine 'EXPLT' to give 
the user one more informal test for expcnentiality . It 
estimates the values of and and if ^ 

and/or ^ 6 then the data may not have an Exponential 
distribution . 

Estimating the values of and , 'EXPLT’ uses 

3 4 

unbiased estimator for (j , S[ (X-^) ] and E( (X-y) ] using 

the formulas: 

n 2 1/2 

a = (.r (X.-X) / (n-1) ) 

1- 1 1 



E[ (X-y) ] = 



n .3 

(X.- X) 
1=1 1 

(n-1) (n-2) 



E[ (X-y) ] 



n _ 4 

(n(n-2)+3) (,^ (X.-X) ) 

1=1 1 



(n-1) (n-2) (n-3) 



3 (7 (n-1) (2n-3^ 

n (n-2) (n-2) 
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2 



structure 



The Subroutine 'EXFLT' is also a fCRTEAN callable 
subroutine with each call producing a plot for the given set 
of data and at the same time calculating and printing the 
estimated values of the parameters "J^and of the data. 

Basically the program is divided into two parts: 

Ihe first part is involved with all required 
computations calculating the Exponential Scores and the 
various statistics needed to compute the ^ and 
parameters. 



Ihe second part is the Subroutine 'EELOT'. ’EEIOT’ 
scales the data points according to their range and plots 
the scaled data along 110 equal-spaced positions of the 
X-axis. No plot is given if the data have constant value or 
if there exists a data point less than zerc. 



A complete description of hew 'EXEIT' opera 
given by the Subroutine and a summary description is 
on the terminal by typing the command DESCBIrS SXPIT 
the ens envirenmenx. When the user types DESCRIBE 
the following response is printed on xhe xerminal: 



tes is 
given 
under 
EXELT 



SDBECUTINE EXPLT 



•EXPIT' takes a s 
increasing erder and plot 
versus exponential scores 
statistics) computed by : 



et of data, sorts 
s the created erder 
(expected values 



it into 
statistics 
of order 
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paramete 

CALL E 

X I 

SCOEES I 
N I 



3. 



plot wi 
that the 
graph i 
feeling 
linear. 
Exponent 
user to 

statisti 
distribu 
for the 
the valu 
6 then 
distribu 
data do 



E[X 



(i) 



] = 



^ ( i/ (n-k+1) ) , 

lc=1 



Also, it computes the estimate 
rs . 

It is called by : 

XELT ( X, SCOEES, N ) 

ARGUMENTS 

s the array containing the data 
s a Mork array of dimension N 
s the number of data values 

Mere information is given in the su 

Interp reti ng The Output 

If data have an Exponential distri 
11 tend to a straight line. But i 
inverse is not generally true. The 
s an indication only and gives 
of the distribution cf data. If t 



bow 
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* EX PLT * with Data Ge aerated fro m V ar i ou s 

Eistr ifcu tioDs 

Ihe Subroutine 'EXPLI’ has been used here, plotting 
data generated by the computer from various distr ibutions, 
in order to get a visual sense of the behavior of the plot. 

Figure 14a gives the plot of ICC Uniform (C,1) 
variates. Obviously nonlinearity governs the plot (as we 
should expect) and thus it suggests departure from the 

Exponential distribution. Besides the nonlinearity of the 
plot the estimates of 'J^and (0.1 and -1.15 res pect ively) 
are far away from the corresponding values (2 and 6) of 
exponentially distributed data. Figure 14b gives the 
associated histogram. 

Figure 15a plots variates generated from a 

Triangular Symmetric distribution (0,2). The nonlinearity 
and the values of the estimates and suggest that one 

reject the exponential assumption. 

2 

Figure 16 shows the plox of X generated variates 

with 10 degrees of freedom. For the same reasons as above 
the exponential model is rejected. 

Figure 17a plots Normally (N (1000,1)) distributed 
generated data and the nonlinearity of the plot is apparent. 
Eesides the nonlinearity, the estimates of and 

(.09 and -0.63) suggest not only to reject the 

Exponentiality of data but in addition lead us to test the 
data for Normality. Figure 17b gives the associated 
histogram. 
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HUflBEB OF OBDEBEO PXIBS » 100 



o 
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lA A 
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Figure 14a - 



A SAMPLE CF SIZE 100 FROM A UNIFORM (0,1) 
DISTRIBUTION 
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Figure 14b - 



HISrOGRAM OF THE DATA Or FIGURE 14a 
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NUMBER Of ORDERED FAIRS = 100 



Figure 15a - A SAMPLE OF SIZE 100 FROM A TRIANGCLAR 
SYMMETRIC (0,2) DISTRIBUTION 
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Figure 15b - 



HISTOGRAM OF THE 



DATA Of FIGURE 
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NUtIBEB or OBDEBEO FAIBS « 100 
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Figure 16a 



2 

X (10 d.f.) BATA 
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Figure 16b - HIsrOGRAM OF THE DATA CF FIGURE 16a 
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HUflBEB OF OBDEBED FAIRS » 100 
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Figure 18a - 



A SAMPLE CF SIZE 50 FROM AN EXPONENTIAL 
DISTRIBUTION 
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Figure 18b - 



HISTOGRAM OF THE DATA Of FIGURE 16a 



75 



IV. SOESgUTINE LIST 



A. DESCBIfTICN 



Another Subroutine presented in this thesis is the 
Subroutine ’LIST'. Its purpose is to list a given set of 
data in ascending order and, taking advantage of like 
occurrences in the data, to print the ordered data in a 
compressed form. This feature becomes highly useful when 



listing a large nu 


mber of data 


points rhat conta 


in 


many 


repeated values. 


It is 


also a tool for 


finding 


multiplicities in 


su pposedly 


continuous data. 


and 


a 



probability function estimating routine for da*a which is 
known to be discrete. 

A complete description of how 'LIST' cperares is given 
in the Subroutine. However a summary is given by typing 



CESCHIBE 


LIST. 


When the user types 


the 


command 


DESCRIBE 


LIST 


under the CMS environment 


the 


following 


response 


will be 


given on the terminal; 







SUBHCOTIKE LIST 

'LIST' sorts a set of data into increasing order and 



gives a 


5-cclumn 


print-out as follows: 




1st 


column: 


Serial number of the 


first occurrence 


this value in th 


e ordered List 




2nd 


column: 


Value of ordered data 


-val ue 


3rd 


column: 


Frequency of occurrence of the value 
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4th cclumn: Percent for the value 

5th ccluan: Graphical representation cf the frequency 

fcr each value. 



It is called by; 

CALI LISI ( X, N ) 
where : 

X Is the array of data 
N Is the numrer of data-values. 

acre information is given in the subroutine. 



INTESPHEIING THE OUTPUT 



The print-cut is a visual representaticn cf the data and 
of each data value frequency. For dara points having the 
same value 'IISI' will print this value once with the number 
cf occurrences. A 5-column output will be printed and its 
interpretation is: 

1. First column gives the serial number of the first 
occurrence of this value in the ordered List 

column gives the value cf ordered data-value 
column gives the frequency cf occurrence of 



2. Second 

3 . T h ir d 
the value 

4. Forth 

5. Fifth 
frequency for each value. 



column gives the percent for the value 
column is a graphical representation cf the 
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E xa mp le : 



Let 1, 1, 


5, 3, 5, 


1, 2, 2, 5, 6 be a given 


data set. 


Then the 


ordered 


data will be: X =X 

(1) 


K) 

II 

X 

II 


X =X =2, 

(^) (5) 


X =3, 

(6) 


X =X =X = 5, X 

(7) (8) (9) (10) 


= 6. Then 


the print-out 


will have 


the follcwing form: 




SEBIAL NOMBEE 


OEDEEED 


DATA FREQUENCY PERCENT 


PHOE. GRAPH 


1 


1 


3 .3 




4 


2 


2 .2 




6 


3 


1 . 1 




7 


5 


3 .3 




10 


6 


1 . 1 




If there are 


no data-p 


oints having a commcn 


value then 


'LIST' gives 


this in 


dication and prinrs cnly 


the ordered 


data. This 


happens 


when data have a 


continuous 


distribution . 








C. USING 'LIST' WITH 


TELEPECNE DATA 1 AND TELEPHONE EATA 2 



Subrcurine 'LIST' was used with Telephone Data 1 
(Figure IS) and Telephone Data 2 (Figure 20) and a brief 
analysis of the output follows. 



Looking at Figure 19 and Figure 20 we may get some 
information of each data set. It can be seen that both of 
the data sets contain a large number of multiple occur lences 
of the data value one and data value two. Is we can see 
reading the fourth column the occurrence of cnes is 1SX for 
Telephone Data 1 and 24% for Telephone Data 2. 
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Comparing also Figures 19a-c, 20a-d we can see that 
multiple cccutrences happen in the range 1 to 24C for 
Telephone Data 1 and in the range 1 to 132 for Telephone 
Data 2. 

Furthermore, guick visual information concerning the 
range where se have multiple occurrences can be obtained 
from the Probability Graph. Thus we can see that for 
Telephone Data 2 there is a region from 113 to 132 where the 
multiple occurrences of values is larger than in neighboring 
regions. Therefore 'LIST' gives the user a useful tool for 
analysis of data. 
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SSRIAL 


QBOE8EO DATA 


PBEQUENCIES 


PERCENT 


PROBABILITT GRAPH 


1 




0, lOOOOOOOE 


01 


128 


0. 19048 




129 




0.20000000E 


01 


54 


0. 080 36 




183 




0. 30000000E 


01 


28 


0.04167 




2 11 




O.aOOOOOOOE 


01 


22 


0.03274 


««*«•«« 


233 




0.50000000E 


0 1 


17 


0.02530 




250 




0.6C000000E 


0 1 


1 1 


0.01637 


« « « 


26 1 




0. 7COOOOOOR 


0 1 


10 


0. 0 1488 


• « « 


271 




0. 80000000E 


0 1 


12 


0.01786 


« * « « 


233 




0. 90000000E 


0 1 


14 


0.02083 


* « « * 


297 




0. 1COOOOOOE 


02 


9 


0.01339 


« « * 


306 




0. 1 1000000E 


02 


10 


0.01488 


« « « 


3 16 




0- 1 2000000E 


02 


1 1 


0.016 37 


« • * 


327 




0. 1 3000000S 


02 


6 


0.00893 


« « 


333 




0, 14000000E 


02 


6 


0.00893 


* • 


339 




C. 15000UC0E 


02 


6 


0. 00893 


« • 


3u5 




0. loCOOOOOF. 


02 


8 


0.01190 


« « « 


353 




0. 17000000E 


02 


8 


0.01190 


« « * 


36 1 




0. 180000COE 


02 


5 


0.00744 


« * 


366 




0. 1 9000000E 


02 


12 


0.01786 


« • « * 


378 




C. 2COOOUOOE 


02 


1 


0. 00 149 




379 




0. 21000000S 


02 


5 


0.00744 


« « 


38U 




C. 2200000UE 


02 


5 


0.00744 


« • 


389 




0. 2 JOOOOUOE 


0 2 


3 


0. 00446 


« 


392 




0. 24000000E 


02 


7 


0.01042 


« • 


399 




0. 250GOG00S 


02 


3 


0. 00446 


« 


U02 




C. GuOOJOOOE 


02 


3 


0.00446 


• 


U05 




0. 27000000E 


02 


2 


0.00298 


« 


u07 




0. 2HOOOOOOE 


02 


3 


0. 00446 


« 


4 10 




0. 29000000S 


02 


5 


0. 00744 


« « 


4 15 




0. 3COOOOOOE 


02 


6 


0.0089 ) 


• * 


421 




0. 3 louooooe 


02 


4 


0.00595 




425 




0. 32000000E 


02 


4 


0. 00595 


« 


429 




0. 3 JOOJOOOE 


02 


2 


0.00298 


« 


4 31 




0. 34000000E 


02 


4 


0.00395 


« 


435 




0. 350J0U1J0E 


02 


3 


0.00446 


« 


4 38 




0. 37000000E 


02 


2 


0.00289 


« 


440 




0. 3M000UUUF 


02 


2 


0.00298 


« 


442 




0. 39000000R 


02 


1 


0.00149 




44 i 




0. UOOOOOODE 


02 


2 


0. 0029H 


• 


445 




0.4 10000005 


02 


2 


0.00298 


« 


447 




0. 4 30000IJOE 


02 


1 


0.00149 




448 




0. 440000005 


02 


1 


0.00149 




449 




0. 450000005 


02 


4 


0. 00595 


0 


453 




0. 4O0000U0E 


02 


3 


0.00446 


« 


456 




0. 47000000E 


02 


1 


0.00149 




457 




0. uaOOOOOOE 


02 


2 


0.00298 


« 


459 




0. 49C00000E 


02 


1 


0.00149 




4o0 




0. 5COOOOOOF 


02 


1 


0.00 149 




46 1 




0, 5 lOOOOOOS 


02 


2 


0. 00298 


« 


463 




0. 520000UOB 


02 


1 


0.00149 




464 




0. 53000000E 


02 


1 


0.00149 




465 




C. 540000005 


02 


3 


0. 00446 




468 




0. 55000000E 


02 


2 


0.00298 


« 


470 




0. 56O000OOE 


02 


2 


0.00299 


« 


472 




0. 5HOOCOOOE 


02 


2 


0.00298 


« 


474 




0. 59000000F 


02 


1 


0.00149 




475 




0. 62000000E 


02 


*) 


0.00298 


« 


477 




C. 63000000E 


02 


3 


0.00446 


« 


480 




0. 64000000E 


02 


1 


0.00 149 




481 




0. 65000000E 


02 


2 


0.00298 


« 


483 




u. Gbooooooe 


02 


1 


0.001 49 




484 




C. 66000000E 


02 


3 


0.00446 


« 


487 




0.69000000E 


02 


1 


0.00149 




468 




0. 7COOOOOOE 


02 


1 


0. 00 1 49 




4S9 




0. 73000000E 


02 


1 


0.00149 




490 




0. 74000000E 


02 


2 


0.00298 


« 


492 




0. 75000000E 


02 


1 


0.00149 




493 




0. 79000000S 


02 


3 


0. 004 46 


« 


496 




0. 33000000S 


02 


1 


0.00149 




497 




0. 84000000E 


02 


1 


0.00 149 




498 




0. 80OOOOOOB 


02 


1 


0.00149 




499 




0 . aaooooooB 


02 


2 


0.00298 


« 


501 




0. 39000000E 


02 


1 


0.001 49 




502 




C. 90000000E 


02 


2 


0.00298 


* 


504 




C. 9 lOOOOOOE 


02 


1 


0.00149 




505 




C. 93000000E 


02 


1 


0.00149 




5 06 




0. 95000000E 


02 


1 


0.00149 






Figure 19a 




USING LIST 


WITH 


TELEPHONE DATA SET 1. 
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507 


0. 98000000E 
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0.00149 


508 


0. 99000000E 
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0.00149 


509 


0. 10600000E 


03 
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0.00149 


510 


0, 1 G900000E 
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511 
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5 15 
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5 16 
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0.00149 


517 


0, 1 1900000E 


03 


1 


0.00149 


518 


0. 12000000S 


03 


4 


0.00595 


522 


0. 12100000E 


03 


2 


0.00298 


524 


0- 12200000E 


03 


2 


0.00298 


526 


0, 12300000H 


03 


4 


0.00595 


530 


0. 12400000S 


03 


3 


0.00446 


533 


0. 1 2800000S 


03 


1 


0.00149 


534 


0. 13500000E 


03 


1 


0.00149 


535 


0. 14200000E 


03 


1 


0.00149 


536 


0, 1 4300000E 


03 


1 


0.00149 


537 


0. 15300000E 


03 


1 


0.00149 


538 


0. 1560000UE 


03 


1 


0.00149 


539 


0. 15800000E 


03 


1 


0.00149 


540 


0, 16100000E 


03 


1 


0.00149 


541 


0. 16500000E 


03 


1 


0.00149 


542 


0. 175000U0E 


03 


1 


0.00149 


543 


0. 17600000E 


03 


1 


0.00149 


544 


0. 1 7700000E 


03 


1 


0.00149 


545 


0, 18300000E 


03 


1 


0.00149 


546 


0. 18600000E 


03 


1 


0.00149 


547 


0. 18700000E 


03 


1 


0.00149 


548 


0. 1920000 OE 


03 


1 


0.00149 


549 


0. 19 300000E 


03 


1 


0.00149 


550 


0. 20200000E 


03 


1 


0.00149 


551 


0, 2 1700000E 


03 


2 


0.00298 


553 


0. 22400000S 


03 


1 


0.00149 


554 


0. 22600000E 


03 


1 


0.00149 


555 


0, 22800000E 


03 


3 


0.00446 


558 


0. 23100000E 


03 


1 


0.00149 


559 


0. 23400000E 


03 


1 


0.00149 


560 


0, 23700000E 


03 


1 


0.00149 


561 


0. 2J900000E 


03 


2 


0.00298 


563 


0. 24000000E 


03 


4 


0.00595 


567 


0. 24 100000E 


03 


1 


0.00149 


568 


0, 24400000E 


03 


1 


0.00 149 


569 


0, 24800000E 


03 


1 


0.00149 


570 


0, 24900000E 


03 


1 


0.00149 


571 


0, 25 lOOOOOE 


03 


1 


0.00149 


572 


0, 25200000E 


03 


1 


0.00149 


573 


0- 27000000S 


03 


1 


0.00149 


574 


0. 27900000S 


03 


1 


0.00149 


575 


0. 28600000E 


03 


1 


0.00 1 49 


576 


C, 29700000E 


03 


1 


0.00149 


577 


0,30300000E 


03 


1 


0.00149 


578 


0, 31200000E 


03 


1 


0.00149 


579 


0- 3 1800000E 


03 


1 


0.00149 


580 


0, 34700000E 


03 


1 


0.00149 


581 


0, 35400000E 


03 


1 


0.00149 


582 


0. 36000000E 


03 


1 


0. 00149 


583 


0, 36400000E 


03 


1 


0.00149 


584 


0. 37000000E 


03 


1 


0.00149 


585 


0,39000000E 


03 


1 


0.00149 


536 


0. 39400000E 


03 


1 


0.00149 


587 


0, 43500000E 


03 


1 


0. 00 149 


588 


0. 46000000E 


03 


1 


0.00149 


589 


0, 46500000E 


03 


1 


0-00149 


590 


0. 47300000E 


03 


1 


0.00149 


591 


0. 48000000S 


03 


1 


0.00 149 


592 


0. 49 lOOOOOE 


03 


1 


0.00149 


593 


0. 54900000E 


03 


1 


0.00149 


594 


0. 6C000000E 


03 


2 


0.00298 


596 


0, 60900000S 


03 


1 


0.00149 


597 


0, 62 lOOOOOE 


03 


1 


0.00149 


598 


0.71100000E 


03 


1 


0.00149 


599 


0. 8 1700000E 


03 


1 


0.00149 


600 


0, 83600000E 


03 


1 


0.00149 


601 


0. 92700000E 


03 


1 


0.00149 


602 


0, 1 1240000E 


04 


1 


0.00149 


603 


0, 1 1500000E 


04 


1 


0.00149 


604 


0, 12700000E 


04 


1 


0.00149 


605 


0, 12890000E 


04 


1 


0.00149 



Figure 19b - USING LIST WITH TELEPHONE DATA SET 1 

(cont . ) 
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606 

607 

608 

609 
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611 

612 

613 

6ia 

615 

616 

617 

6 18 

6 19 

620 

621 

622 

623 

624 

625 

626 

627 
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631 

632 

633 

634 

635 

636 

637 

6 38 

639 

640 

64 1 

642 
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658 
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664 
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668 

669 
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19c 



0. 12980000E 


04 
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0.00149 


0, 13050000E 


04 
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0.00149 
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04 
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0.00149 


0, 134800002 
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0.00149 
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0.00149 
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0.00 149 
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0.00149 


0, 15190000E 


04 


1 


0.00149 


0. 15470000S 
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0.00149 


0. 16330000E 


04 
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0.00149 


0. 17870000S 


04 


1 


0.00149 


0. 20720000S 


04 


1 


0.00149 


0. 24830000E 


04 


1 


0.00149 


C. 23060000H 


04 


1 


0.00149 


0. 29620000E 


04 


1 


0.00149 


0. 30260000E 


04 


1 


0.00149 


0. 3^0 10000E 


04 


1 


0.00149 


0. 3593000UE 


04 


1 


0.00 149 


0. 3b850000S 


04 


1 


0.00149 


0, 395200GOE 


04 


1 


0.00149 


0. 4 1570000S 


04 


1 


0.00149 


0. 44690000E 


04 


1 


0.00149 


0. b20dOOOOE 


04 


1 


0.00149 


0. 76 140000E 


04 


1 


0.00149 


0. B3220000E 


04 


1 


0.00149 


0. 9C 15000CS 


04 


1 


0.00149 


0. 96250000E 


04 


1 


0.00149 


0. 98060000E 


04 


1 


0.00149 


0. 981BOOOOE 


04 


1 


0.00149 


0. 10154000E 


05 


1 


0.00149 


0, 1039BOOOE 


05 


1 


0.00149 


0. 1045100 GE 


05 


1 


0.00149 


0. 1 0939000E 


05 


1 


0.00149 


0, 1 1280000S 


05 


1 


0.00149 


0, 1 34470 0 0E 


05 


1 


0.00149 


0. 14385000E 


05 


1 


0.00149 


0. 15155000E 


05 


1 


0.00149 


0. 15294000E 


05 


1 


0.00149 


0. 15504000S 


05 


1 


0.00149 


0, 1 5847000E 


05 


1 


0.00149 


0. 15868000S 


05 


1 


0.00149 


0. 16280000E 


05 


1 


0.00149 


0. 16299000E 


05 


1 


0.00149 


0. lb361000S 


05 


1 


0.00 149 


0, 1 640300 0E 


05 


1 


0.00149 


0. 16817000E 


05 


1 


0.00149 


0. 17 174000E 


05 


1 


0.00149 


0, 17667000E 


05 


1 


0.00149 


0, 18218000E 


05 


1 


0.00 149 


0. 18649000E 


05 


1 


0.00149 


0. 19461000E 


05 


1 


0.00149 


0, 2 1848000E 


05 


1 


0.00149 


0. 23493000S 


05 


1 


0.00149 


0, 24692000E 


05 


1 


0.00149 


0. 26443000E 


05 


1 


0.00149 


0. 3C974000E 


05 


1 


0.00149 


0. 35644000E 


05 


1 


0. 00 1 49 


0, 38003000E 


05 


1 


0.00149 


0. 40131000E 


05 


1 


0.00149 


0. 47120000E 


05 


1 


0.00149 


0. 47592000E 


05 


1 


0.00149 


0, 6 1710000S 


05 


1 


0.00149 


0. 69775000S 


05 


1 


0.00149 


0. 85993000S 


05 


1 


0.00149 
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SEBIAL SaSBEa 



0BC8B20 DATA 



FHEQ0EMCI2S 



PEBCSIfT 



PB03ABILITT GBAPH 



1 


0. 10000000E 


01 


178 


0. 24185 




179 


0. 20000000B 


0 1 


36 


0. 0489 1 




215 


0.30000000E 


0 1 


n 


0.01495 


• • 


226 


0.4G000O00B 


01 


6 


0.00815 


• 


2S2 


0.50000000B 


01 


6 


0.00815 


• 


238 


0. 60000000B 


01 


5 


0.00679 


« 


243 


0.70000000B 


0 1 


5 


0.00679 


* 


248 


0. 80000000B 


01 


4 


0.00543 


« 


252 


0. 90000000B 


01 


4 


0. 00543 


« 


256 


0, lOOOOOOOR 


02 


9 


0.01223 


«« 


265 


0. 1 1000000B 


02 


2 


0.00272 




267 


0. 12000000B 


02 


3 


0.00408 


* 


270 


0. 13000000B 


02 


1 


0.00 136 




27 1 


0. 14000000B 


02 


1 


0.001 36 




272 


0. 15000000B 


02 


1 


0.00136 




273 


0. 1»3000000B 


02 


4 


0.00543 


* 


277 


0. 19000000B 


02 


1 


0.00136 




278 


0.21000000B 


02 


1 


0.00 1 36 




279 


0. 220000002 


02 


1 


0.00 1 36 




280 


0.24000000E 


02 


3 


0. 00408 


• 


283 


0. 2500000CE 


02 


3 


0.00408 


• 


286 


C. 26000000B 


02 


2 


0.00272 




238 


0. 2700000UE 


02 


1 


0.00136 




289 


0. 3COOOOOOB 


02 


3 


0.00408 




292 


0. 32000000E 


02 


1 


0.001 36 




29J 


0. 36000000S 


02 


1 


0.00 1 36 




294 


0. -lOUOOOOOB 


02 


2 


0.00272 




296 


0. U2000000S 


02 


1 


0.00 1 36 




297 


0. 44000U00B 


02 


1 


0.00136 




298 


0. 470000002 


02 


4 


0.00543 


* 


302 


0. 480000 JOB 


02 


1 


0.00136 




303 


0.49000000S 


02 


9 


0.01223 


• « 


312 


0. 500000J0B 


02 


4 


0.00543 


« 


316 


0. 53000000R 


02 


3 


0.00408 


« 


319 


0.55000000B 


02 


2 


0.00272 




321 


0. 560000002 


02 


4 


0.00543 


* 


325 


0. 57000000B 


02 


3 


0.00408 


m 


328 


0. 5a000000B 


02 


1 


0.00136 




329 


0. 59000000B 


02 


3 


0.00408 


* 


332 


0.6GOOOOOOE 


02 


2 


0.00272 




334 


0. 62000000E 


02 


3 


0.00408 


* 


337 


0.63000000B 


02 


1 


0.001 36 




338 


C. 6U000000E 


02 


a 


0.0054 3 


« 


34^ 


0. oSOOOOOOS 


02 


3 


0.00408 


m 


345 


0. 66000000B 


02 


Q 


0.00543 


m 


349 


0,670000008 


02 


1 


0.00136 




350 


0. 6HOOCOOOB 


02 


5 


0.00679 


m 


355 


0. 690000002 


02 


3 


0.00408 


m 


358 


0. 70000000B 


02 


3 


0.00408 


• 


361 


0. 7 10CD000B 


02 


5 


0.00679 


« 


366 


C. 73000000E 


02 


1 


0.00136 




367 


0. 74000000E 


02 


1 


0.00136 




368 


0.75000000E 


02 


1 


0.00136 




369 


0. 76000000B 


02 


1 


0.001 36 




370 


0. 77000000S 


02 


1 


0.001 36 




37 1 


0. 81000000S 


02 • 


1 


0.00 1 36 




372 


0. 86000000E 


02 


1 


0.00136 




373 


0. 880000002 


02 


1 


0.001 36 




374 


0. 89000000E 


02 


1 


0.001 36 




375 


0. 93000000B 


02 


1 


0.00136 




376 


0. 95000000B 


02 


1 


0.001 36 




377 


0, 96000000B 


02 


1 


0.00 1 36 




379 


0. 97000000B 


02 


1 


0.00 1 36 




379 


0. 990000O0B 


02 


3 


0.00408 


• 


382 


^ 0.10000000B 


03 


1 


0.00136 




383 


0. 10 lOOOOOB 


03 


1 


0. 001 36 




384 


0. 10500000B 


03 


1 


0.00136 




335 


0, 10600000E 


03 


1 


0,00 1 36 . 




336 


0. 10700000B 


03 


1 


0.0013b 




387 


0. 1C900000E 


03 


1 


0.00136 




338 


0. 1 10000 008 


03 


1 


0.00 1 36 




389 


0. 1 1 lOOOOOB 


03 


1 


0.001 36 




390 


0. 1 1 200000B 


03 


1 


0.00136 




391 


0, 1 1300000E 


03 


4 


0. 00543 


0 


395 


0. 1 1U00000E 


03 


3 


0.00408 


« 


398 


0. 1 1500000B 


03 


1 


0.00 1 36 




399 


0. 1 1600000E 


03 


4 


0,00543 


« 


‘ 4 03 


0. 1 1700000B 


03 


1 


0. 00136 
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404 

409 

4 17 

426 

433 

445 

453 

460 

469 

475 

478 

480 

484 

487 

490 

492 

493 

494 

495 

496 

497 

498 

500 

501 

503 

504 

505 

507 

508 

509 

5 10 

511 

512 

513 

514 

515 

516 

517 

5 18 

519 

520 

523 

525 

526 

528 

530 

532 

533 

534 

536 

538 

539 

540 

541 

542 

543 

544 

545 

546 

547 

548 

549 

550 

552 

553 

554 

555 

556 

557 

558 

559 

560 

56 1 

562 

564 

565 

566 

567 

568 

569 



0. 1 1300000E 


03 


5 


0,00679 


♦ 


0. 1 1900000B 


03 


8 


0.01087 


** 


0. 12000000B 


03 


9 


0,01223 




0. 12100000R 


03 


7 


0. 00951 




0. 12200000E 


03 


12 


0.01630 




0. 12J00000E 


03 


8 


0.01087 




0. 12400000E 


03 


7 


0.00951 




0. 12500000E 


03 


9 


0.01223 




0. 12600000E 


03 


6 


0.00815 


*■ 


0. 12700000E 


03 , 


3 


0.00408 


* 


0. 12B00000E 


03 " 


2 


0.00272 




0. 12900000E 


03 


4 


0,00543 




0. 13000000E 


03 


3 


0.00408 




0. 13200000E 


03 


3 


0.00403 


* 


0. 13300000E 


03 


2 


0,00272 




0. 13700000E 


03 


1 


0.00136 




0. 14000000E 


03 


1 


0.00136 




0. 14300 00 OE 


03 


1 


0,001 36 




0- 1520000 OE 


03 


1 


0.001 36 




0. 15800000E 


03 


1 


0.00136 




0. 16600000E 


03 


1 


0,00136 




0. 16900000E 


03 


2 


0.00272 




0. 17000000E 


03 


1 


0.001 36 




0. 1 730000 OE 


03 


2 


0.00272 




0. 17600000E 


03 


1 


0.00136 




0. 17800000S 


03 


1 


0.00136 




0. 18000000E 


03 


2 


0.00272 




0. 13200000E 


03 


1 


0.001 36 




0. 18500000E 


03 


1 


0.00136 




0. 18700000E 


03 


1 


0.00136 




0. 19000000E 


03 


1 


0.00 136 




0. 19400000E 


03 


1 


0.00136 




0. 19500000E 


03 


1 


0.00136 




0. 19900000E 


03 


1 


0,00136 




0. 20600000E 


03 


1 


0.00 1 36 




0. 20900000B 


03 


1 


0.00136 




0.2 1600000E 


03 


1 


0,00136 




0. 22900000E 


03 


1 


0.00136 




C. 23000000S 


03 


1 


0.00136 




0. 23700000E 


03 


1 


0.00136 




0. 2 J900000E 


03 


3 


0.00408 




0. 24000000E 


03 


2 


0.00272 




0. 24100000E 


03 


1 


0.001 36 




0. 24400000E 


03 


2 


0.00272 




0. 24700000E 


03 


2 


0,00272 




0. 24800000E 


03 


2 


0,00272 




0. 25000000S 


03 


1 


0.00136 




0. 25100000E 


03 


1 


0.00136 




0. 25400000E 


03 


2 


0.00272 




0. 25500000E 


03 


2 


0.00272 




0. 25600000E 


03 


1 


0.00136 




0. 25800000E 


03 


1 


0.00136 




0. 27200000E 


03 


1 


0.00136 




0. 27500000E 


03 


1 


0.00136 




0. 28000000S 


03 


1 


0,001 36 




0. 29800000E 


03 


1 


0,001 36 




0. 3 I800000E 


03 


1 


0,001 36 




0. 32400000E 


03 


1 


0,001 36 




0. 34000000E 


03 


1 


0.001 16 




0. 34400000E 


03 


1 


0.00136 




0. 34600000E 


03 


1 


0.00 1 36 




0. 35900000E 


03 


1 


0. 00136 




0. 36500000E 


03 


2 


0,00272 




0, 36700000E 


03 


1 


0.00 1 36 




0. 37700000E 


03 


1 


0.00136 




0. 38000000E 


03 


1 


0.00136 




0. 38700000E 


03 


1 


0.00136 




0. 4 1300000E 


03 


1 


0,00136 




0. 4 1800000E 


03 


1 


0.00136 




0.42900000E 


03 


1 


0. 00136 




0. 45300000S 


03 


1 


0.00136 




0. 4o600000E 


03 


1 


0.00136 




0. 47200000E 


03 


1 


0,00136 




0. 50300000E 


03 


2 


0.00272 




0. 5C400000E 


03 


1 


0,001 36 




0. 50700000E 


03 


1 


0.00136 




0. 5 lOOOOOOE 


03 


1 


0,00136 




0. 51400000E 


03 


1 


0.001 36 




0.51700000E 


03 


1 


0.00 1 36 




0. 52800000E 


03 


1 


0.00 136 
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570 


0, 57000000E 


03 


1 


0.00136 


571 


0,57200000E 


03 


1 


0.001 36 


572 


0. 58600000E 


03 


1 


0. 001 36 


573 


0. 58700000E 


03 


1 


0.00 1 36 


574 


0. 60100000E 


03 


1 


0.00136 


575 


0. 6 1400000E 


03 


1 


0.00136 


576 


0. 62400000E 


03 


1 


0.00136 


577 


0, 64000000E 


03 


1 


0.001 36 


578 


0. 66300000E 


03 


1 


0.00136 


579 


0. 71400000E 


03 


1 


0.00 1 36 


580 


0, 71600000E 


03 


2 


0.00272 


582 


0. 74000000E 


03 


1 


0.001 36 


583 


0. 74200000B 


03 


2 


0.00272 


585 


0. 79 100000B 


03 


1 


0.00136 


586 


0. 85300000E 


03 


1 


0.00136 


587 


0. 86200000E 


03 


1 


0.00136 


588 


C. 86600000E 


03 


1 


0.00136 


589 


0. 87500000E 


03 


1 


0.00 1 36 


590 


0, 38500000E 


03 


1 


0.00136 


591 


0. 89400000E 


03 


1 


0.001 36 


592 


0. 92600000E 


03 


1 


0.00136 


593 


0. 92800000E 


03 


1 


0.00136 


594 


0, 94100000E 


03 


1 


0.00136 


595 


C. 976000U0E 


03 


1 


0.00136 


596 


0. 99100000E 


03 


1 


0.00136 


597 


0. 1021000 0E 


04 


1 


0.00 1 36 


598 


0. 10230000S 


04 


1 


0. 00136 


599 


0, 1C620000E 


04 


1 


0. 00 1 36 


600 


0, 1U860000E 


04 


2 


0.00272 


602 


0, 1 1070000E 


04 


1 


0.001 36 


603 


0, 1 1 180000E 


04 


1 


0.00136 


604 


0. 1 1 25000 OE 


04 


1 


0.00136 


605 


0. 123 10000E 


04 


1 


0.00136 


606 


0, 1 2500000E 


04 


1 


0. 00 1 36 


607 


0. 12610000E 


04 


2 


0.00272 


609 


0. 12960000B 


04 


1 


0.001 36 


6 10 


0. 13670000E 


04 


1 


0.001 36 


611 


0. 14 120000E 


04 


1 


0.00136 


612 


0. 14130000E 


04 


1 


0.001 36 


6 13 


0. 14620000E 


04 


1 


0.00136 


6 14 


0. 14850000E 


04 


1 


0.00136 


615 


0. 14970000E 


04 


1 


0.00136 


616 


0. 1 5040000E 


04 


1 


0.00136 


6 17 


0. 15170000E 


04 


1 


0.00136 


6 18 


0. 15J40000E 


04 


1 


0.00136 


6 19 


0. 1 535000 CE 


04 


2 


0.00272 


621 


0. 15620000E 


04 


1 


0.00136 


6 22 


0. 15940000E 


04 


1 


0.00136 


623 


0. 15950000E 


04 


1 


0.00136 


624 


0. 16030000E 


04 


1 


0.00136 


625 


0. 16520000E 


04 


1 


0.00136 


626 


0. 165J0000E 


04 


1 


0.00136 


627 


0. 16930000E 


04 


1 


0.00136 


628 


0. 17000000E 


04 


1 


0.00136 


629 


0. 17150000E 


04 


1 


0.00136 


630 


0. 17500000E 


04 


1 


0.00 1 36 


631 


0. 17550000E 


04 


1 


0.00136 


632 


0. 17860000E 


04 


1 


0.00136 


633 


0. 1B220000E 


04 


1 


0-00136 


634 


0. 1824000 OE 


04 


1 


0-00 1 36 


635 


0. 18550000E 


04 


1 


0. 001 36 


636 


0. 18770000E 


04 


1 


0.001 36 


637 


0. 18930000E 


04 


1 


0.00136 


638 


0. 18970000E 


04 


1 


0.00136 


639 


0. 19230000E 


04 


1 


0. 001 36 


640 


0, 1 9640000E 


04 


1 


0.00136 


641 


0, 19650000E 


04 


1 


0.001 36 


642 


0. 205 10000E 


04 


1 


0.00 1 36 


643- 


0. 20630000E 


04 


1 


0.00136 


644 


0. 21670000E 


04 


1 


0.00136 


645 


0. 21850000E 


04 


2 


0.00272 


647 


0. 22080000E 


04 


1 


0.001 36 


648 


0. 22530000E 


04 


1 


0. 00136 


649 


0, 22b00000S 


04 


1 


0.001 36 


650 


0. 22710000E 


04 


1 


0.00136 


651 


0, 22870000E 


04 


1 


0. 00 1 36 


652 


0. 231 60000E 


04 


1 


0.00 1 36 


653 


0. 23210000E 


04 


1 


0.001 36 


654 


0. 24290000E 


04 


1 


0.00136 


655 


0, 24390000E 


04 


1 


0.00 1 36 
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656 


0* 24680000E 


04 


1 


0,00136 


657 


0. 24720000E 


04 


1 


0,00136 


658 


0, 249000002 


04 


1 


0,001 36 


659 


0, 26840000E 


04 


1 


0.00136 


660 


0, 27930000E 


04 


1 


0,00136 


661 


0. 28820000E 


04 


1 


0,00136 


662 


0, 289100002 


04 


1 


0.00136 


663 


0. 29380000E 


04 


1 


0.001 36 


664 


0. 30800000E 


04 


1 


0,00136 


665 


0. 3 1 170000S 


04 


1 


0.00136 


666 


0, 316500002 


04 


1 


0.00 1 36 


667 


0. 32860000E 


04 


1 


0.00136 


668 


0. 33150000E 


04 


1 


0.00136 


669 


0, 337800002 


04 


1 


0.00136 


670 


0. 35100000E 


04 


1 


0.00136 


671 


0. 38070000E 


04 


1 


0.001 36 


672 


0. 3037000 OE 


04 


1 


0.00136 


673 


0,3d500000E 


04 


1 


0.001 36 


674 


0, 38580000E 


04 


1 


0.00136 


675 


0. 38790000E 


04 


1 


0.00136 


676 


0, 38880000E 


04 


1 


0.00136 


677 


0. 394400002 


04 


1 


0.00136 


678 


0.397800C0E 


04 


1 


0,00136 


679 


0. 43730000S 


04 


1 


0.001 36 


680 


0, 43930000E 


04 


1 


0.001 36 


681 


0,451200002 


04 


1 


0.00136 


682 


0. 46370000E 


04 


1 


0.001 36 


683 


0. 475900002 


04 


1 


0.00136 


684 


0, 4835000 CE 


04 


1 


0.00136 


685 


0.52160000E 


04 


1 


0.00136 


686 


0, 54920000E 


04 


1 


0.00136 


687 


0, 56520000E 


04 


1 


0. 00136 


688 


0. 57540000E 


04 


1 


0.00136 


689 


0. 6 1350000E 


04 


1 


0.00136 


690 


0. 624 10000E 


04 


1 


0.00 1 36 


691 


0. 63720000E 


04 


1 


0,00136 


692 


0,63850000E 


04 


1 


0.00136 


693 


0. 66590000E 


04 


1 


0.00136 


694 


0. 68160000E 


04 


1 


0.00136 


695 


0. 60210000E 


04 


1 


0.00136 


6 96 


0. 73070000B 


04 


1 


0.00136 


697 


0. 73290000E 


04 


1 


0.00136 


698 


0. 77460000S 


04 


1 


0.00136 


699 


0. 79270000E 


04 


1 


0.00136 


700 


0, 8039000 OE 


04 


1 


0.00136 


701 


0, 80580000E 


04 


1 


0.00136 


702 


0. 83530000E 


04 


1 


0.00136 


703 


0. 88470000E 


04 


1 


0,00136 


704 


0. 92060000S 


04 


1 


0.00136 


705 


0, 92560000B 


04 


1 


0.00136 


706 


0, 95 170000E 


04 


1 


0.001 36 


707 


0. 954 10000E 


04 


1 


0.001 36 


708 


C. 95620000E 


04 


1 


0.00136 


709 


0. 9965000 CE 


04 


1 


0.00 1 36 


7 10 


0, 1C020000E 


05 


1 


0.001 36 


711 


0. 10376000E 


05 


1 


0.00136 


7 12 


0. 1C837000E 


05 


1 


0.001 36 


713 


0, 10918000E 


05 


1 


0.00136 


7 14 


0, 12042000E 


05 


1 


0,00136 


715 


0, 12793000S 


05 


1 


0.00136 


716 


0. 13012000E 


05 


1 


0.00136 


717 


0. 1 31 79000E 


05 


1 


0,00136 


718 


0. 14692000E 


05 


1 


0.001 36 


719 


0. 15825000E 


05 


1 


0.00136 


720 


0, 16877000E 


05 


1 


0.00 1 36 


721 


0. 18885000E 


05 


1 


0, 00 136 


722 


0. 19273000E 


05 


1 


0.00 136 


723 


0. 19696000E 


05 


1 


0.00136 


724 


0. 19849000E 


05 


1 


0.00 1 36 


725 


0. 1985 lOOOE 


05 


1 


0.00136 


726 


0, 2C846000E 


05 


1 


0,00136 


727 


0, 2 1440000S 


05 


1 


0.00136 


728 


0. 24573000E 


05 


1 


0.00136 


729 


0. 247700002 


05 


1 


0.00136 


730 


0. 26278000E 


05 


1 


0.00136 


731 


0. 27238000E 


05 


1 


0.001 36 


732 


0. 281 J3000E 


05 


1 


0.001 36 


733 


0. 29913000S 


05 


1 


0.00 1 36 


734 


0. 39867000E 


05 


1 


0.00 1 36 


735 


0. 53138000E 


05 


1 


0. 00 1 36 


7 36 


0, 6727100 0E 


05 


1 


0.00 1 36 



Figure 20d - USING LIST WITH TELEPHCNE DATA SET 2 

(cont . ) 
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V. 



ASSE SSME NT OF VAEIAB I LIT Y BOURNES 



A. miECDOCTION 



In th€ previous sections Prcbatility Plotting 



Subroutines were presented for use in informal estimation of 
the form of the distribution of a set of independent 
observations. In the present section subroutines computing 
the variability of some basic statistics are described. 

Me know that J. (sample mean) is an unbiased estimator of 

the population mean ( = x f (x) dx ) , where X is a 

X 

population with unknown distribution (? (x) ) and unknown 



observations from X. Furthermore we oan compute the 
Variance of 1 as: 



X 



mean, and X , X 



2 



/ • • • f 



X is a sample of independent 
n 



X. 




1 




n 



n 



(J 



2 



Var[ X] 



n 



n 
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Thus W6 see it is possible to estimate the Var[X], even 

2 

if the population Variance ( C ) is uckncwn, by simply 

2 2 

usinvj the sample ^variance S as follows; Var[I] = S /n, 
where S = Z (X -X) / (n-1) (an unbiased estimator of CJ~) . 






Therefore S /n is an unbiased estimator of Var[I] also. 



Furthermore, especially for a Normal population, we can 
obtain a confidence interval for the mean using the 
t-statistic . 



Eut, although for I there exists a direot assessment of 
variability, for estimates of other population parameters 
such as Skewness, Kurtcsis, Coefficient of Variation, and so 
forth, this is not so simply obtained. Thus several methods 
have been introduced to obtain assessments of variability 
and for two of them (Sectioning cf data and the Jacknife) 
the Subroutines 'SECTN' and 'JACK' will be presented. 



E. SUBfiCaTISE SECTN 



1 . Eescri pt ion 

'SECTN' Subroutine is used for assessing variability 
of estimates from data based on the Sectioning Method. 

I 
A 

. . . , X 

n 

function 

8(n) = 9 



he basic idea of this method is: 

ssume we have n independent observations X , X , 

1 2 

from a population with unknown distribution 

F (X). Let & be a parameter cf F (x) and 

X . . , . , -X o 

(X ,X ,...,X ) be a statistic which estimates u . 
1 2 n 
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Now W€ want to estimate the variance of 9(n) which is a 
new Random Variable and we work as follows: 

1. Divide the sample into r disjoint sections, of 

size m (r should be such that n=mr. If this is not possible 
some of the latter data is discarded) . ^ 

2. For each section form the same estimate (m) , 

i~1#2/.««/r 

3. Cempute the average of 9{C'^),that is; 

- 1 r ^ 

d (n) = ^ 9 (m) , 

r 1 = 1 1 

doing so, we have an estimate of 9 and 

- 1 r 

E[9(n) 3 = E[ .X 9. (m) ] 

r 1=1 1 

1 

= — rE[9(m) 3 
r 

/VJ 

= E[ 9(m) 3 , 

also we have: 

— 1 IT 
Var[9(n) 3 = Var[-— Z 8 (m) ] 

r i= 1 i 

r 

Var[ 2T 9 (m) ] 
i=1 i 

Var[ 9 (m) 3 

1 r — 2 

Z (d. (m)-9(n)) 

r-1 1=1 1 

2 
S- 

9(m) 
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Ihus we see that 7ar[ 9 (ni) ] and Var[ 9 (n) ] can be 

unbiasedly estimated from the sample variance of Bi(ai) , the 
main advantage and the purpose of the sectioning of data 
method. 



Eur the main disadvantage of rhis method is that we 

would like r to be as large as possible in order tc make 

2 

the variability of the variance esrimate S as small as 
possible. This however may be worse for the bias 
properties of the estimation procedure. Therefore the 
choice of r is a factor that should be considered. 



It should be noted here that if (m) are 
approximately normal variates (This can be tested by the 
previous described Subroutine 'NORMPL') then confidence 
intervals for the unknown parameter can be obtained based 
on the t-Statistic, in the following way: 



5 ^ 9(n) 

y{n) ± t 

1/2 (1-^/2) , (r-1) 

r 

where ; 

9 (n) is the mean of the sectioned data statistics 
obtained from the column named 'mean' of the second table 
('estimated parameters of the sample parameters') of the 
program output. 

1/2 

is the standard deviation of the 



^ 9(n) 

sectioned data statistics divided by the sguars root of the 
number of sections, obtained from the last column of the 
same table of the program output. 



(1- ac/2) , (r-1) 



is the (1- <^/2) quantile of the 



t-distributicn with r-1 degrees of freedom. 
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2 . ££ 231^5 S tr uctur e 



*S2CIN*, using the Sectioning method, estimates the 
following statistics: Mean, Median, Variance, Standard 
CeviaticD, Coefficient of Variation, Skewness, Kurtosis, 
Minimum, and Maximum. (The formulas which have been used to 
compute these estiaates are described in the comments of the 
Subroutine.) The first table is then printed bj the 
program, containing the values of these parameters. Then 
•SECTN*, using the computed estimates of all sections for 
each parameter, estimates the Mean, Median, Variance, 
Skewness, Kcrtosis, and Standard Deviation divided hj the 
sguare root of the number of sections. The second table is 
printed containing these values. 

There are three restrictions in using 'SECTN': 

1. The number of sections must be no greater than 100. 
If it is, a diagnostic message is printed and only estimates 
from unsecticned data will be given. 

2. The number of data values must be greater than 3; 
otherwise a diagnostic message is printed without any 
calculation. 

3. The size of each section must be greater than 3. If 
it is less than, or egual tc 3, then the program gives 
estimates fcr the entire set of (unsectioned) data. A 
diagnostic is printed. 

A ccmplete description of how 'SECTN' operates is 
given in the Subroutine. Furthermore, a Summary is given by 
typing on the terminal the command DESCRIBE SECTN under the 
CMS environment. The following response will be given on 
the terminal when the user types the above command: 



91 



SOBHCUTINE SECTN 



•SECTN* is intended tc estimate a set of tasic 
statistics cf a given set cf observations using the 
•Sectioning cf data method'. Also for each estimated 
statistic, estimates of some basic statistics such as the 
mean, standard Deviation, and so forth, are given. 

It is called by : 

CALL SECIN ( X, N, K ) 

where : 



X 


Is 


the 


arra 


y of 


data 








N 


Is 


the 


na 


mbe r 


cf 


data values 


(must 


be greater 




th 


an 3) 














K 


Is 


the 


numb 


er of 


desi 


red sections 


(no 


great er 



than 100) 



Note: k should be a number which minimizes the number of 
data points that will have to be discarded. 'SECTN* places 
the data into the equal size sections discarding any data 
left over. 



Fcr k<3 or k>1C0 or 
unsectioned data will be 
estimated statistics will be 



(n/k)<3 only estimates from 
given and no estimates fcr the 
computed . 



Nc output is expected if n<3 . 

here information is given in the subroutine. 
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•3 • Usina * SEC TN I with Tel e ph one pa^ Set 2 

'SBCIil' was used on Telephone Data Set 1 to assess 
the variability in the Mean, Median, Variance, Standard 
Deviation, Coefficient of Variation, Skewness and Kurtcsis. 

The 672 data-points of Telephone Data Set 1 were 
broken down into 16 disjoint sections with 42 data-points 
per section. Because of this break-down no data-pcints were 
discarded. Compare the values of row 'unsectioned' (see 
figure 21) with the values of the corresponding statistics 
computed by the Subroutine ' HISTGS/HISTFS ' (see figure 1) . 
The values are the same. 



Now if we want to assess the variability in ary of 
the above parameters we proceed as follows: 

1. We take the mean (0 ) of the parameter, whose the 



variability 


ue 


wanr 


to assess, from the 


tab 


' ESTIMATED 


EABAMETEHS 


OF SAMPLE PA5AMETEES'. 




2. 


We ra 


ke. 


from 


the same table, the value 


(S) 


last column 


(STD. 


DEV. 


/NS**. 5) for the same paramet 


3. 


Using 


the 


t-sta 


tistic with k-1 d.f. and 


the 



9 ± St 

(1-^2), (k-1) 

we get a (1-<?~)% confidence interval for rhe paramerer. 



Example 

In order to find a 95% confidence inrerval for the 

Skewness of Telephone Eaxa 1 we take: 

8= 4.9343, S = .30312, t = 2.131 

. 975, 15 



Therefore the 95% confidence interval is: 

4.9343 ± .30312*2.131 = 4.9343 ± .64595 
===> [ 4.2884, 5.58C2]. 
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It should be acted here that the use of the Variance 
estimate from the sectioned data in order to get confidence 
intervals of the parameters is based on the normality 
assumption and the independence of the estimates from the 
sections. The normality will depend on the number of 
data-points in each section, which should be kept large. 
This reguirement , however, conflicts with the need to make 
the number of sections large to reduce the variability in 
the estimate of the variance of the statistics. The 
skewness estimates of each of the 16 sections, can be run 
through the normal plotting routine to see whether the use 
of the t-statistic confidence interval is valid. 
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Figure 21 - USING SECIN WITH TELSPHCSE DATA SET 1. 
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SKElrfNESS 4.9343E 00 4.7465E 00 1.4701E 00 1.2125E 00 2.4573E-01 -1.0259E-01 -1.4647E 00 3.0312E-OI 

KURTOSSIS 2.45S2E Ol 2.19bOE 01 1.2484E 02 1.117JE 0 1 4.550HE-OI 4.0534E-03 -K5801E 00 2.7933E 00 



c. 



SUBROUTINE JACK 



1 • Descri ptio n 

This routine is used for assessing variability of 
estimates of data based on the Jacknife method, anc for 
reducing the bias in estimates also. It is particularly 
useful fcr data with small sample size. 



A big picture cf the method has as fellows; 

Ler I , X , ..., X be a sample of n independent, 
12 n 

identically distributed observations from a population with 



unknown distribution function P (x) . Alsc let S te an 

X 

unknown parameter of P (x) to be estimated. Further mere 

X 

suppose a method (biased or unbiased) is available for 



estimating 0 . Then we proceed as follows; 

1. Divide the sample size inro r disjoint groups, each of 
size m (r should be such that n = mr. If this is not 
possible seme of the last data-points will be discarded.) 

2. Compute the estimator 9 (n) of 9 based on all n = mr 

ebservatiens . ^ 

3. Compute the estimator 9 (n-m) based on the n-m 

i 



th 

ebservarions , having deleted the i group. 

A. 

4. Compute the so called 'PSEUDO VALUES' d (n) ; 

i 

9 (n) = r9(n) - (r-l)9 (n-m) , i = 1,2,. ..,r . 

i i 
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Then the Jacknife estimator is defined to te the average of 
Eseudo Values: 



A 1 r 

9(n) = 2 : 8 (n) 

r 1 = 1 1 



1 2 ^ IT ^ 

( r 8(n) - (r-1) S’ 8 (u-m) ) 

r i= 1 i 



= r 8(n) - ----- z 9 (n-m) 

r i= 1 i 



Havirg so defined the Jacknife estimator it can be 

proved that this estimator is an unbiased estimator cf 9, 

-2 

except for terms n higher order, assuming that the 

bias has tne form cf E[ 9(n) ] = ^ + a/n + terms of higher 

crder in n. That is, the Jacknife estimator eliminates a 
-1 

n bias term. Namely if 

E[ 8 (n) ] = 9 + an + 0 (n ) , 

then we have: 

A j-_1 j- ^ 

E[8(n) ] = E[r8(n) Z 9 (n-m) ] 

r i= 1 i 



rE[0(n) ] .X £,[ 9 . (n-m) ] 

r i=1 i 



o -1 -2 r-1 (\ ra r -2 

r ( ^ + an + 0 (n ) ) (r>7+ + C (n )) 



n-ii n-m 



r8+ ---- + ...-(r-1) (8+ +...) 

n m (r- 1 ) 



= r 



Q ra a (r-1) a Q 

y + + , + '^ +. 

n m (r-1) 
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I 




mr m 



= 9 . 



A special case 
mentioasd here, called 
Khere r=n and in=1. 



of Jacknife estiaator can be 
the 'Complete Jacknife Estimator*, 



Seme properties of the Jacknife esticator follcw: 

-2 

1. It is unbiased up to order n 

2. The 'rseudo values* can be used to obtain variance 

estimates for the Jacknife estimator since they can be 

considered (see £18]) as approximately independent and 

-1 r 2 

identically distributed. Thus (r(r-1)) ^ (8 (n)-0(n)) 

i=1 i 

A. 

should be an approximate estimate of Var[9 (n) ] , 
and 



(9(n) -S) (r(r-l)) 



1/2 



t a Q 21/2 

£. 2 . (0 (n) - y(n)) } 

1 = 1 1 



should be approximat elly distributed as a t-statistic with 
r-1 d.f. This procedure is particularly useful if the 
number n cf data-points is small, but it must be used with 
care . 

3. For large n it can be shown that for the ccaplete 
Jacknife estimator, under very general conditions we have: 

A 

Var[ 8 (n) ] > Var[ 8 (n) ] . 



98 



■■’li 



m 



.S 




2. 112311™ Structure 



Subroutine 'JACK* is a FOETRaN-callable Suhrcutine 
which takes a set of data, groups it into r disjoint groups 

and for each set of the r-i group computes and prints in a 
table the statistics; Mean, Median, Variance, Standard 
Deviation, Coefficient of Variation, Skewness and Kurtosis. 
The same statistics are also computed for the ungroupec set 
of data. Then, using these statistics, the Jacknife 

estimator of the above parameters are computed and printed 
in another table along with its Variance and Standard 
Deviation . 

Ihe program is divided into the main program ard tne 
Subroutine JACKES. The main program groups the data and 
successively calls the Subroutine 'JACKES' in order to get 
the estimates of the above parameters for ungrouped data and 
for each group as well. A table is then printed containing 
the estimated parameters. The 'Pseudo-Values ' for each 
parameter are also computed by the main program, and are 
used to compute the Jacknife estimator for each parameter 
and its Variance and Standard Deviation. 

Tc avoid division by zero the number of observations 
must be greater than three, otherwise no output is given. 
For the same reason the expression (r-1) (n/r) must evaluate 
tc greater than three, otherwise the program will give only 
estimates of ungrouped data. 

A complete description of how 'JACK' operates is 
given in the Subroutine. A Summary is printed on the 
terminal by typing DESCRIBE JACK under the CMS environment, 
which responds; 
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SUBHCUTINE JACK 



Jack is intended to estimate the statistics : Mean, 
Median, Variance, Standard Deviation, Coefficient of 
Variation, Skewness and Kurtosis of a given set of 
independent observations, using the Jacknife method. In 
addition to the above parameters which are given for each 
group, the Jacknife estimator along with its variance and 
standard deviation is given for each parameter. 



It is called by: 



CALI JACK 
where : 

X 

XS 

3TAT 

N 

IG 



( X, XS, STAT, N, IG ) 

is the array of data of dimension N 

is a work array of dimension N (returns ordered 

data) 

is a two-dimensional work array, dimensioned by 
(IG,7) 

is the number of data-values 

is the number of groups plus one (=r+1) 



Note : r must be such a number as to rainimiz 
number of data-pcints that will have to be discarded 
places the data into the equal size groups discardin 
data left over (last observations in order of or 
input) . 



€ the 
. JACK 
g any 
iginal 



No output is expected if N < 3. Furthermore if 
IG < 2 or (IG-2) (N/(IG-1)) < 3 only estimates for ungrouped 
data will be given . 

More information is given in the Subroutine. 
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3. 



SsifiS * JAC K * w ith Teleph one Da ta 1. 

To assess the variability in the Jlean, Median, 
Variance, Standard Deviation, Coefficient of Variation, 
Skewness and Kurtosis of Telepfone Data 1 the Sabrcutine 
•JACK* is DOW used. In order to avoid to discarding any of 
the data-foints we use 16 groups with 42 data-points per 
group . 



Observing figure 22 we see that the values of the 
parameters for ungrouped data are exactly the same as the 
corresponding values of the same parameters of figure 1 
which was froduced by using the Subroutine HISTGS/HISTFS. 
Now in order to assess the variability in any of the above 
parameters we have to use the t-statistic with (r-1) 
d.f. along with the values printed under 

•ESTIMATED JACKNIFE PASAMETEHS' using: 

9 ± St 

, (r-1) 

where 9” is the Jacknife estimator and S is its Standard 

Deviation. For example to assess the variability of the 

Skewness with a confidence level .05 we have: 

9= 7.37321, S = .900118, t = 2.131 , 

.975,15 

Therefore a S5% confidence interval of the Skewness is: 
7.37321 ± .900118^2.131 = 7.37321 ± 1.S1815 

===> [5.4551, 9.2914] . 

This should be compared with the point estimate 4.9343 and 
the confidence interval estimate [4.2884, 5.5802] obtained 
from subroutine SECTN. The data is so skewed that one will 
inevitably have trouble here with the Jacknife procedure. 
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Figure 22 - USING JACK WITH TELEFHCNF DATA SET 1 
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KUBTOSSIS 6.70793E 01 14.6U056B 03 1.71037E 01 



VI. 



singijj^ebvee Fi8si::cgt3E fissi;Seev|d ggsuEs «ith 

»EAR Mfi» S TRUCTU RE 



A. THE EABEA{P,Q) FRCCESS 



The starting point for these processes is the defirition 
of a first-crder autoregressive model for a stationary 
sequence of random variables {X _ } ; 



X — rX i — 

i i- 1 i 



If the marginal distribution of the X is fixed to be 

i 

exponential with parameter Tv for all i; 



-JVx 

P {X < X} = 1-e , Tv> 0, x> 0 , 

i 



then 6. is zero with probability r and exponential (Ov) with 



pronability T^r. Thus 



X 

i 



rX 


w. p. 


r 


i-1 






•rX + E 


w. p. 


1-r, 


i-1 i 







( 1 ) 



where {E_} is a sequence of i.i.d. exponential (2\) random 
variables. The process defined by (1) is the exponential 
autoregressive process of order 1 the EAR(1) process. The 
ccrrelaticn structure is r(j) = r , j = 0,1,2,... . 
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The first order moving average exponential process, 
EilA(1) is defined as 

' 'Si. , V .p. S 

(0<^<1;i = 0,±1 ,±2, . . .) 

I + E w.p. ^^-S) 

i i- 1 



and has 


correlations ^(1) = 


0-S) 


and p(1) = 0, j = 2,3,... 


The 


general EMA (g) model 


takes 


the form 





w.p. b 


q 1 


q+1 


^ i E 


w.p. b 


g 1 g-1 x-1 


q 



/e + / E +...+£ E 

i ^q-1 i-1 1 i-g+1 



( 2 ) 



w.p. 



b 

2 









£ 

'9 i 



+ 




E 

i-1 



+ 



+ 



S E + 

1 i-g+1 



E 

i-g 



w.p. 



h 



1 



for i = 0, ±2, 

1 2 g 

A 

b. = ^ (1- ^) . . . (1-^. ) 4 

1 g X 1 - 1 

( 1 - if ) ... (1 ) 

q i 



where 

i = q+1 f 
q ^ i ^ 2, 
i = 1 . 



The serial correlations are given by the equation 



(q) 

F 



(j) 



corr (X ,X ) 
i i-r 



If 



^ n b 

V= j V V + j 



< 



0 



(1 ^ j ^ q) 



(g+1 < j <Qo) 
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Tc ccQvert this EMA (q) process into a mixed 
autoregressive moving-average exponential process of orders 



1 and g. 


called EMA (1 , g) we 


replace E 


in 


(2) with 








i 


-q 






rA 




w. p . 


r 




A = •< 


i-q-1 










i-q 


rA 


+ E 


w. p. 


(1-r) 






i-q- 1 


i-q 




For 


further 


results 


on these 


processes and their 



properties see Gaver and Lewis (1978), Lawrance and lewis 
(1977), Jacobs and Lewis (1977), and Lawrance and Lewis 
( 1978) . 



E. USE CE AUTOREGRESSIVE (EAR(P)), MOVING AVERAGE (EKA(q)) 
AND MIXEE STRUCTURES ( EARM A ( P , Q) ) IN MOEELLING QUEUES 



Consider for simplicity a queue with a single input 

stream and a single server, and a f irst-ccme-f irst-served 

(EIFO) service discipline. Let S , i = 0,1,2,... , denote 

i 

t h 

the service time for the i arrival, and let X , 

i 

th 

i = 1,2,... , denote times between arrival of the i and 



(i-1) customers. As is usual we assume that the first 



customer (with service time S ) arrives at time zero and 

0 

finds the queue empty. 



If the {S } and {X } sequences are i.i.d. exponential 
i i 

random variables with parameters ^ and o. respectively, we 



have the M/M/1 queue. 
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L j! 



iyHiii 



' - ! 












Now let 

E be exponential (;^) and independent, i = 0, + 1,... ; 

i 

g be exponential (a) and independent, i = 0, + 1, . 

i 

We want to model queues with correlated (autocorr elated 
and/or cross-correlated) service and inter-arrival times, 
the service and inter-arrival times both having marginally 
exponential distributions. We also want the gueing model to 
include the a/M/l queue as a special case. 

There are five simple possibilities which we pur forward 
here, based on the use of EARHA(p,q) processes, giving what 
we call an EAfiMA aD/MD/I queue. 

In what follows let p and p be, respectively, the 

S X 

order of the a ut or egressive components of the service time 

sequence {S } and the inter-arrival time sequence {X } , and 
i i 

let q ,q be, respectively, the order of the moving average 
S X 

component of {S } and {X } . These parameters can take 

i i 

values C,1,... . If p = 0, q = 0, the sequence is 

independent; if p = 0, the process is purely moving average 
and if q = 0, then the process is purely autcregresive . 

The five possibilities or cases are as follows: 

1) Let IS } be EARaA(p ,q ) over {£ , S , ...}; 
i S S i i-1 

let IX.} = {£.}. 

1 1 

Thus the arrivals are a Poisson process and the service 
times are autocorr elated , i.e. a mixed autoregressive 
moving average exponential sequence. , 



106 



2 ) 



1 



Let S = E, i = C,1,2, ; 

i i 

let X te EABMA(p ,g ) ever {6 f S > — }• 
i XX i i-1 



Then the service times are independent and the arrival 
process is non-Poisson because of the dependency between 
the inter-arrival times, i.e. the arrival process is a 
point process with EARMA(p,g) structures. 



3) Let {S } 
i 


te EARMA (p , g ) 
S S 


ever 


1 


2. r 

i- 1 


. . . J . 


Let {X } 
i 


te EARMA (p , g ) 
X X 


ever 


ie.. 

1 


^i-r 


. . . ^ . 


The serv 


ice and arrival 


processes 


are 


autocorrela ted. 


but the 


two processes 


are 


inde 


penden 


t. The marginal 


distribu 


tions of {S } an 
i 


d {X.} 
1 


are 


still 


exponential . 



4) To ccuple the two processes, with resultant dependence 
in the arrival and service processes, the sinplest 
procedure seems to be the following. 

Let {S } be EARflA(p ,g ) in the following sense: 
i S S 

{S } is EARflA(F ,g ) over {E., ^ • 

i S S 1 > i > 1-1 

Then if X = & , i = 0,1,2,... we have that {S } is an 
i i i 

autcccrr elated seguence and also cross-ccrrelated with 

{X = 6 }, i.e. (S , X } is a bivariate dependent 
i i i i 

seguence cf random variables with expenential marginal 

distributions, although {X } is still Poisson. Note 

i 

that the S seguence is an autccor related seguence 
i 

because of the cross-coupling, but not a pure 
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It does have an exponential 



EARMA (p ,q ) sequence. 

S S 

marginal distribution hovever. 



As an example let p = 1, q = 1. 

S S 

Let 



S = E 
0 0 



S = ^ E 
1 S 1 


w .p. 




A = 
1 


CX 


6 

1 




= ^ E + A 
S 1 1 


w ,p . 


(1- 


s 








s = ^ E 
2 S 2 


w .p* 




A = rA 
2 1 






w . p . 


= / E + A 
2 2 


w .p . 


(1-- 


^ ) = rA 

S 1 


+ 


2 


w. p. 



S 

i 



/ E 



w . p . 



w . p . 




r A 

i-1 




w. p. r 


rA 


+ 


w.p, (1 


i-1 


i 



In these equations one can replace 6 with X and then 

i i 

the cross-coupling between the sequences {S } and {X } 

i i 

becomes apparent. 



Inter pretatio n . We have positive correlation between S 

i 

and q previous inter-arrival times (when p = 0) . If 
S S 
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i-p -1) 

s 



the S s (j = 

j 



are short, thee S 

i 



will he short if all the q previous inter-arrival times 

S 

are short. This models the c ase where the ser ver t end s 
M£ g ets lon^. Cf cource he also 

slows down when the queue gets shert and it is not 
immediately clear what the effect on an average waiting 
time will he. 



4‘) Similar to 4 but {X } is EABHA(p ,g ) over 

i XX 

^ e 2 l2 , 2L? /••• ands = E. 

i cx i OL i-1 i i 

Interpretation is net clear, but it would have the 

same effect as balking in the input stream. If service 

times get long (and presumably the queue gets long) then 

inter-arrival times get long. 

5) A mere general possibility which allows one to model 
dependenoy in the input stream, in the service stream 
and dependency which couples them is as fellows: 

Let {D } be a unit exponential independent sequence, a 
i 

driving sequence. 



Let 


1^.} 

1 


be 


EASMA(p ,q } 
S S 


over {E , E ,. 

i i- 1 




Let 


1^.} 

1 


be 


EARMA (p , g } 
X X 


ever {a./ &. /• 

1 1- 1 


..} . 


Let 


{S.} 

1 


be 


EAEMACp ,q 

C^ C^ 


) over {E , , ID . , 
1 1 


1- 1 


L 


{X.} 

X 


be 


EAEMA<p ,q 


) over (I , o-D , 
i i 


f m • • ^ • 

i- 1 
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This is a very general scheme which reduces tc all or 
the previous schemes as special cases. Fcr example if 



F =g =p =g =0, then we have case 3. If 

Cl Cl Cj_ 

q =1/ F = 0, / = 0, P = 9 = 0 ; p = 0 (and p 

C^_ Cj Cj, S S X Cl 

and g are not equal to zero) we have case 4. 

In all cases the basic equation for M/H/1 queue waiting 

time H still holds: 
n 

M =0 
0 



S =(W+£-X ) n=1,2,.... 

n+1 n n n+1 

These equations can be used to generate successive W s 

n 

in a simulation. We have not touched on hew the correlated 
sequences are started. This is somewhat arbitrary as it is 
not known (except for the EMA(1) and MA(1) processes) hew to 
start them so as tc produce a stationary exponential 
sequence. The problem is aggravated, as in cases 4, 4 and 

5, when we are dealing with bivariate exponential sequences. 
Then the aarginal processes may be stationary but not the 
bivariate process. 
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C. PROGEAM STEDCTURE 



1. General 



Ic simulate the EABMA(p,g) model, a BORTRAN Program 
has been written to take care of all the existing cases of 
that model and also of the a/M/1 queue. Cepending on the 
input values which the parameters (as described below) will 
take, the program can be used for creating in each run a 
particular case of the EARMA(p,g) queue. Eecause of that 
generality cf the program it is not suggested that it be 
used for these EARMA cases where either the arrival tines or 
the service times cr both are not autccorr elat ed or 
cross-correlated. The reason for this suggestion is 
efficiency. It is suggested for use with the EARM£(p,q) 
cases where both arrival and service times are 
aurocorrelated and, furthermore, cross-correlated with any 
order. The program is divided into the main progran and 
into the Subroutines BETAS, AOTOR and EAREA. 



2* Msin Pr ogram 

The main program performs all I/C operations and 
calculations required for the desired statistics. 
Specifically it is designed: 
a) To read 

1. All randem number generator seeds required to 
generate the sequences of the Exponential and Uniform 
variates to be used during rhe execution. 

2. The values for the variables; 

N : Eumber of arrivals to be generated. 



Ill 



(1 fcr the 



a : Sumter of replications. 

RX, RS : Arrival rate, service rate. 

CIX, CFS, CPG : Taking the values 1 cr 0 

coupled process) . 

KXl, FX, KS1, KS : The order of mcviEg average for 

arrival and service processes. 
Taking the values 1 or 0 (1 for 

autoregressive case) . ^ 

Arrays with the values of for 

moving average part. 

RH0X1, RHOX, HHCSI , RHCS : Parameter values for the 

autor egressive parts. 



QXl, QX, QSl, Q£ 
EX, BX1, BS, BS1 



t) To calculate and print 
SUflX (I , J) ; 

sots (I, J) ; 
w (l,J) ; 

WE (I,J) ; 

D (l,J) ; 

fcr I = 1,2, 3,4, 
where : 

(In what follows k gets the values N/4, N/2, 3N/4, N 
fcr 1=1,2, 3,4 respectively and J= 1 , 2 , . . . , M) 



SUtlXM (I ,J) ; 
SOMSM (I,J) ; 
wa (I, J) ; 

WME (I, J) ; 

EM (I,J) ; 

J = 1,2,.. 



• t ^ t 



k (J) 

SGMX (I, J) = Z X , 
r=1 r 



X 



(J) 



the r 



th 



arrival for 



th 

the J replication in tne BARBA gueue; 

k (J) (J) th 

SUEXM(I,J) = z (XM) , (XB) the r arrival 

r= 1 r r 

th 

for the J replication in M/M/1 Queue; 



k (J) (J) th 

SOBS (I, J) = ZS , S the r service for 

r=1 r r 

th 

the J replication in the EARMA gueue; 
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J. (‘J) (J) th 

SOMSil{I,J) - ^ (St3) , (SM) the r service 

r= 1 r r 

th 

for the J replication in M/M/1 Queue; 

(«J) (J) (J) (J) 

a (I,J) = = max{ (W + S - X ) , 0} 

k k-1 k-1 k 

th th 

= viaxtxng txme of k arrival in J replication 
of the correlated queue; 

(•J) 

WM|I,J) - (WM) The same as W(I,J) but 

k 

for M/M/1 queue; 

k (J) 

WB (I, J) = 27 w /k 

r=1 r 

X (J) 

WM£{I,J) = X (WM) /k 
r=1 r 

D (I,J) = {SUMS (I, J) - SUMX(I,J)}/K - 1/RS + 1/5X 

DM|I,J) = (S0MSM(I,J) - SUMXM (I, J) }/K - 1/fiS +1/RX 
k (J) 

SDMSM(I,J) = X (SM) , 
r= 1 r 

Mete: Between the gro up- variables (KX1, KX, KSI, KS) , (CXI, 

QX, QS1, QS) , (CPX, CPS, CPG) , (RHOX1, HKCX, EHOSI, HHCS) 
and (BX1, EX, BSl, 3S) there are the following relations: 

1. The number of elements of each variable (array) from 
E-gioup should be the same as the value of the corresponding 
variable frci K-group. 

2. The 'O' value of a variable from Q-greup implies that 
the corresponding variable from H-group is not needed. 

3. The value '1' of CPG dominates any value of CPX, or 
CPS. 

Any valid combination in the values of the above 
variables specially of the Q-group, C-group and K-greup (0 
or not 0), generates a particular EARMA(p,g) case and is 
accepted by the program. 
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The basic approach in the main program is tc 
generate at once and for each replication three independent 
exponential sequences with parameters XL (arrival rate), SL 
(service rate) and 1 (unit exponential) which are stored in 
the arrays EXPXL, SXFSI and EXP1 respectively. Ther the 
Subroutine AOTOR which is involved with the autoregressive 
part or the model is called (if it is reeded for the 
particular EABMA case) . A Icop follows which is executed as 
many times as the value of N. From the loop the EAEMA 
Subroutine is called which is involved with the iicving 
average part of the model. The number of calls of Subrcutine 
EARMA depends on the particular case which is being rur. The 
statistics and all infcrmaticn that will help us to analyze 
the model are also computed inside the Iccp and stored in 
the arrays described above. The program continues execution 
until the desired statistics for all d replications have 
been calculated and gathered. Then the output part of the 
program fcllcws which gives us the values of calculated 
statistics on paper and on punched cards fcr analysis. 

3. Subrcutine BETAS 



return an 
fellows; 



This is a simple Subrcutine and its purpose is tc 
array SUMBXI/X/SI/S with element- values as 



SUMBX (I) = X EX(r) , 
r=1 



I = 1 , 2 , . . . , value of K-greup 



to be used for the choice of the order in the moving 
average. It is called at most 4 times during the execution 
and only if the K-group variables have value greater than 

•0 • . 
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4 



• Subroutine A OTCR 



•AOTCB' accepts in each call an exponential secuence 
(generated previously) of variates and transforms it into an 
autocorr elated sequence. It is called only if any of the 
Q-group variables equals 1. For each Q-group variable having 
value 1 it is called i! times. 

5* BAfi MA 



This is the main Subroutine of the program and it 
has been created so that it can be used for all EASaA(p,q) 
cases. Eecause of the generality quite a fev parameters are 
transferred and in each call just one value is returned. 



This value is the 



th th 

I arrival or the I service time 



as 



it 



has been modified 



called 


only if 


the 


moving 


average ; 


that 


is grea 


ter than 


0. 



because of the loving averag 
particular EABMA(p,g) case 
is only if any of the K-group v 



e. It is 
requires 
ariables 



6 . 



3 is® require ments 



Seme 


CPU 


-times have been gathered i 


n the c 


our se 


the simulation 


of some EABMA(p,q) 


case 


s, usin 


g the 


360/67, FCRIEAN 


H compiler as 


follows : 








N 


M 


K-G30UP QS1 


C-GBOOP 




TIME 




2000 


500 


0 1 


0 


9 


Hin. 51 


Sec . 


5000 


500 


0 1 


0 


22 


Min. 25 


Sec . 


10000 


500 


0 1 


0 


43 


Min. 30 


Sec . 


1000C 


500 


0 1 


CPS=1 


55 


Min. 20 


Sec . 



of 

IBM 
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from these results we see that the value of 11 is the main 
factor that affects the CPO-time (linearly) . We see also 
that one of the C-group variables may increase the CPC-time 
by 255?. We ignore the factor of M (it is another main 
factor as N is) because we may stop the program at any 
replication, getting the last random number generator seeds 
and then contxnuing another time using these seeds. 



D. ANAIXSIS OF SIMULATION FESOLTS 



Since nc analytical properties of the EAHMA(p,q) models 
can be derived, a simulation has been done to srudy their 
properties. But as we have noticed, this program is a 
general program accomodating all EAfiMA(p,g) cases. Since 
our study is limited to only the two cases below, in order 
to save CPU-time a modified program has been usee that 
simulates these particular two cases. 



1 . fiueue wi th D epend e nt Se r vice 



The icdel; 



{X.} = 6. 

1 X 

(S } is EARMA (1 , 0) 

i 



over { E , E , . . . } 
i i- 1 



That is: 



r A 



i-1 



. rA + E 
i-1 i 



w . p. r 



w.p. (1-r) 



where A = E 
0 0 



To simulate this model we chose the values .25, .50, 
.95, .99 for traffic intensiry t and the values .25, .50, 
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.90, .95, .96 for the correlation r, in order to cover a 
representative range. In what follows the string S#tt#rr 
stands for the run with t=tt and r=rr . That is the run 
S#25#50 stands for the run with t=. 25 and r=.50 For each t 
and r we ran the program and the reguired statistics were 
gathered, analyzed and plotted using the Subroutines 
HISTGS/FS, NCRdPL and EXPLT. The sample statistics from 
HISTGS/FS were also tabulated separately frcm the figures so 
that one could obtain an overall picture of whether or not 
the waiting times had converged in mean and in distribution 
to the limiting distribution. 



From the analysis of all runs and 
following results: 

1. Simulating the type S#25#rr and Si5C 
r =.25, .50 and .90) it was possible fc 

reach the steady state for M=2000. On the 
the type SiJ25#rr and S#50#rr (with rr = . 
necessary tc go up to N=10000, in order for 
reach the steady state. Thus we see 
ccrrelaticn affects the choice cf the value 
N to increase as the correlation is inc 
cnly the correlation affects the choice of 
intensity affects it much more, since for 
and S#99tfrr the W and do not reach the s 
N=10000 even if r is low. Because 
requirements we restricted ourselves tc a d 
the particular case S#99#98 increasing 

value of N up to 320,000, where the W and W 

close tc the steady state. Thus we may c 
traffic intensity and/cr high correlation 
value of N in order to achieve steady st 
through 23d, 24a-24c, 25 and 26a-26d ju 

ccnclusicns. In figures 23a-23d, where we 



plots we ha 



#rr model 
r the M and 
other han 



95, .98) 
the W and 
that the 
cf N, req 
reasing. B 
N. The t 
the type S 
teady stat 
cf the C? 
etailed stu 
successive 



appear 

cnclude tha 
require a 
ate. Figur 
stify the 
are dealin 



ve the 

(with 
TI to 
c for 
it was 
to 
high 
uiring 
rt not 
raf fic 
#95#rr 
6 for 
0-time 
cy of 
ly the 

to be 

t high 
large 
€s 23a 
above 
g with 
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th 

the waiting time of the 2000 (W ) arrival for the case 

2000 



t=0 . 50 , r=0 . 25 , 



we can see that W has converged to a 

2000 



value of about 6.3, and furthermore a straight line appears 

in the plot under EXPLT. We also see approximate 

convergence for the correlated queue for the TJ , which is 

2000 

presented cn figures 28a-28c, and its value has converged 
to a value of about 3.03 (figure 23d) . Note that the 



r=0.98) . 
carry th 
stationary gusue. 



ion 


of 


correlation 


into 




the s 


er vice 


tim 


€ 


has 


t h€ 


average waiting 


ti me 


fr 


offl 2.5 


tc a 


bo ut 


3 


.03, 


roxi 


mately 20%. 
















€ ca 


nno t 


say the same 


fer 


figures 


24a- 


24 c 


w h 


ere. 


cf 


t he 


high ccrrelaticn 


(r 


=0.98) 


we 


do no 


t 


have 


ce 0 


f W i 


and H, although N= 


10 


,000, 


and for 


fi 


gur e 


bee 


ause 


cf the high 


traf f 


ic 


inte 


nsity 


(t=0. 


S9 


) we 


V€ c 


cnvergence of W a 


nd ¥. 




Figur 


es 26 


a-26d 




also 


the 


no 


n-cen vergence 


of 


w. 


H, because 


of th 


e 


high 


cn 


and 


the high 


traffic 


int 


ensi ty 


(t 


= 0 


.99, 


Th 


is 


is the case 


for w 


hich we 


eventually 


had zo 


sia 


ulat. 


icn out to about 


N = 


500, 00 


0 to 


obse 


rv 


e a 
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Figure 23a - QUEUE WITH 
SEQUENCE AND POISSON INPUT. 



WITHOUT 



ZEROS 



FROM 



EAR1 AUTOREGRESSIVE SERVICE TIME 
HISTOGRAM OF THE WAITING TIMES 
THE RUN S#50#25; m=500 



2000 

REPLICATIONS, RX=.2, RS=. 4. No OF ZEROS = 255. 
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Figure 23t - QUEUE WITH EAE1 AUTOHEGEESSIVE SERVICE IIHE 



SEQUENCE AND POISSON INPUT. EXPONENTIAL PICT (EXPLT) CF THE 



WAITING TIMES W WITHOUT ZEROS 

2000 

111 = 500 REPLICATIONS, RX = .2, BS = .4. 



FROM THE RUN S#50#25; 
No OF 2EECS=255. 
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figure 24a - QUEUE WITH EAE1 AUTOREGRESSIVE SERVICE 



SEQUENCE ANE POISSON INPUT. HISTOGRAM OF THE WAITING 



W WITHOUT ZEROS FROM 

10000 

REPLICATIONS, RX=2.5, RS= 1 0 ; 



THE 



RUN 



S#25#98 ; 



No OF ZEROS=354. 



TIME 

TIMES 

m=500 
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Figure 24t - QUEUE 
SEQUENCE AND POISSON 
iiAITING TIHES H 

1 00 

01 = 500 REPLICATIONS, 



WITH EAR1 AUTOREGRESSIVE SERVICE TIME 
INPUT. EXPONENTIAL PLCT (EXPLT) CF THE 
WITHOUT ZEROS FROM TEE RUN S#25#98; 

0 0 

RX=2.5, RS=10; No OF ZEBOS=354. 
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Figure 24c - QUEUE WITH EAH1 AUTOREGRESSIVE SERVICE TIMS 

SEQUENCE ANC POISSON INPUT. NORMAL PLOT (NORMPL) CF THE 

WAITING TIMES W WITHOUT ZEROS FROM THE RUN S#25#98; 

10 0 00 

01 = 500 REPLICATIONS, RX = 2.5, RS=10; No OF ZEROS=354. 
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figure 24d - QUEUE Hlia EARl AUTCB EGRESSIVE SERVICE llKES AND PCISSCH INPUT. TAbUlATICN OF SAP.IIE STATISTICS fOR THE 
CISTFItUIICNS Of CUMULATED INTERARRIVAL TIMES AT N = 5000 AND N=* 10000 (X 2 AND X 4), CUMULATED SERVICE TIMES (S 2 AND S 4J , 
WAITING TIMES, kITU CE kllHOUI ZERCS (HI ♦ 0, HI, I FOB CASE N=250C, 5000, 7500, 10C00) , CUMULATED kAITING TIMES (HI BAH etc.) 
AND THE AVERAGED DIFf£RE^CES BETWEEN S2 AND X2 (D2); OEIAJWED FHCM id = 500 BEELICAIICRS CF THE RON Si25«SE; BX=2.5; RS = 10 . 
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Figure 24e - QUEUE 
AND POISSON INPUT. 
WAITING TIKES ¥ 



1 0 0 00 

REPLICATIONS, RX=2.5, ES=10 



WITH EAR1 AUTOREGRESSIVE SERVICE TIMES 
HISTOGRAM OF THE CUMULATED AND AVERAGED 
FROM THE RUN S#25#98; m=500 
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figure 24f - QUEUE WITH EAE1 AUTOHEGHESSI VE SE3VICE TIMES 



AND POISSON 


INPUT. EXPONENTIAL 


PLOT 


(EXPLT) 


OF 


THE 
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AVERAGED WAITING TIMES 


TI 


FROM 


THE 


RUN 
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S#25#98; ai=500 
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Figure 24g - 
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AND AVE3AGEE WAITING TIMES 



QUEUE WITH EAB1 A UTOR EGBESS I V E SERVICE TIMES 
INPUT. NOEMAI PICT (NORMPL) CF THE CUMULATED 

FROM THE 
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m = 500 REPLICATIONS , RX=2.5, RS=10 . 



RUN S#25#98; 
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figure 25 - QOEUE WITH EAB1 AOTOfcEGfiESSIVE SEDVICE HUES AND PCISSCN INfUI. lADUlATICN Of SAfElE STATISTICS fOR TEE 
ClSIflfOTICNS Of CUEULAIED INTEBAflRIVAl TlttES AT N*5000 AND N= 10000 (X 2 AND X U), CUEULAIED SERVICE XlttES (S 2 AND S 4), 
WAITING IIBES, kITH CE kJXHOUT ZERCS (WI ♦ 0, WI, I FOB CASE N=2500, 5000, 7500, 10CC0) , CUMULATED WAITING TIKES (K1 BAB etc.) 
AND TBE AVERAGED DlffEBEECES BETWEEN S2 AND X2 (D2) ; OBTAINED FRCM 0=500 REIIICATICNS CF THE RUN SA99I25; RX=2.97; RS=3 . 
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Figure 26a - 
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Figure 26b 
AND POISSON INPUT 
TIMES W 



WITHOUT ZEROS 

1 000 0 

REPLICATIONS, RX=2.97; RS=3 



QUEUE WITH EAR1 AUTOREGRESSIVE SERVICE TIMES 
EXPONENTIAL PLOT (EXPLI) OF THE WAITING 
FROM THE RUN S#99#98; m=500 



No OF ZER0S=19. 
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QUEUE WITH EAR1 A UTORSG HESS I V E SERVICE TIMES 
INPUT. NORMAL PLOT (NORMPL) OF THE WAITING 
WITHOUT ZEROS FROM THE RUM S#99#93; m=500 

RX=2.97; HS=3 ; No OF ZER0S=19. 



133 



CCKRELATED CUEUE 





cr* 


in <d 






o 






1 














sf 


3 


O 


X 




o 


1 






1 






1 






n o 






o 




















»3 


O' 


X 


X 


n 


ri 




















n rvj 


V 


03 


3 




















O 


O 


X 


n 


X 


3 








1 








r 


<Vi 


-r -Vi 


n 


A 


3 






1 












s#- 


3 


3 


3 


o 


j- 


3 












I 






fO 


fM 


M 




O 




















O' 


o 


O* 


O' 


•n 


o 








i 




I 




'X 






• « 


« 


• 








J 












X 






« 
























r\ 


7 


M 












j 














3 




fM 


_4 


n 


n 


















f\i 




• n 






o 






1 














iM 




m 






O' 


















<Ni 




! M 




































si 
















wO 


o 




o 


*n 


c 


















03 


r- 


m 




O' 


iC 


3 
















o 




>1 


o 


y' 




o 


















o 


fM 




\n 


3 


O' 


O 




















tn si' 






o 




















r- 


Ui 


in 


3 


u 


O 












1 




•* 


nT 


>o 


o 


vO 


l-v 


















*■ 


>0 


fii 


sr 


3 


si 


3 




















c 


'T 


fM 




3 




















x 


O 


s/* 


O' 


3 


O 
















>r 


• 






t 


























• 




• 
























S' 


rf> 






o 


















1 


3 




<M 




vn 


O 
















X 






O 






o 


















1 X 


IM 




fM 






o 












1 






IM 




fM 






♦n 






























m 












' 






CJ 


in o 




n 


o 






i 












j 


n 


r- 


vn 




3 


o 


















rsj 


>» 






O' 


o 






1 














O' 


iil> 


|43 


X 


IM 


3 


j 
















«J 


(• 


'T 


ni 


'T 


o 




















r- 


O' 


4" 


fl 


o 


O 


t 














m 


U’ 






m 


(M 


o 






■ 














3 




wd 


>r 


IX 


o 


1 


















r* 


O 


fM 


03 


f-> 






I 












1 


fM 


O' 


M- 


3 


a; 


o 


j 




1 


1 








1 


• 






• 


• 




































1 




























1 












1 


fM 


ti 


O 




—4 


3 


1 
















i/' 






























fM 




»M 






O' 


















-* 










-r 






1 
























si 
















wO 




fNJ 


»n 


r- 


nT 


o 




m 


V 


n 


o 


f'- 


o 




1 

03 


O 


3 


Vn 


3 


o 


o 


1 


O' 


ct> 


1 X 




X 


o 


o 








fM 




o 






n 


:▼? 






O 




I o 


3 


o 


fM 


X 


3 


3 




J* 


3 


, <n 


X 


si 


o[ 




<»> 


m iM 


o 




c 




in 


n 




—4 


—4 


o 






O 


>r 


03 


r- 


o 


O 




si 






o 




o 




'j- 


^1 o 


•o 


o 


o 


>f 


03 


o- 


lA 




O 


o 






X 




«»> 


sf 




o 


si 




C\t 


rd 


vn 


rd 


o 




O 


■c 


in 


fM 


03 


o 








fM 




X 


o 






O' 


O' 


sT 


3 


X 


o 




3 


o 


si 




X 


Ol 


f<\ 


• 






• 


• 


• 




o 


o 


r*'' 


n 


fM 


o 




-fl 


• 


• 


* 


• 


4 


• 


c 


3 


O 


o 




fs- 


ot 




CO 










o 




o 


o 


o 


fM 


-4 


o 








o 


O 


_d 


_d 


o 




3 


o 


o 


IM 


.4 


31 


X 


03 










o 








• 




t 






X 


rM 




fM 






o 












4 






— • 




-* 






iTW 




o 


o 


O 


o 


O 


o 
















vn 




G 


G 


G 


o 


o 


Gi 


















t 










o 




















1 










ol 




O 




fM 


o 


-iJ 














tn 










sT 


f- 


d4 


o 






1 








vnt 








»n 




i^> 


3 




















O' 


fM 


h.. 


X 


X 


3 






j 










o 


o u 


<r 




o 


















( 




fia 


X 




y3 


o 














1 


r\l 


'T 


IH- 


fM 


n 


■U 


o 


















I ^ 


fM 


X 


O 


sT 


r«. 


o 














i 




r\ 




O' 


fM 


T 


o 




















—t 




n 


sT 


3 


G 
















J, 




4 




• 


• 


• 




h- 


o 


O 


r- 


o 


o 




[ X 


• 


• 


• 




4 


• 




si 


X' 


_ 


o 


X 


o 








m 












<o 


't 


h3 


•o 


o 






O' 


o 


X 




fM 


fM 




fs» 


si 


fM 




X 


o* 




«rv 




fO 












03 


n 


—X 




o 






—4 




<-d 






n 


1 


li' 


O'. 


• in 


n 


si 


oi 








—4 






r< 


IM 


03 


'J' 


>r 


O 


o 


o 
















si 


<M 




fM 


n 


si 


n 


o. 


















fM 






0 




3 


















1 


3 


3 


X 


O' 


ds 


3 


yO 


cr\ 


<M LT 


'J* 


'T 


o 


a 


O 


o 




Q 


fM 


CJ 




03 


-d 


X 


dd 


U 


fVJ 


o 


c 


3 


O 


G 


X 


Zd 


O: 


O 


m 


<r 


m 


03 


—4 


o 




o 


o 


Ci 






f-i 




o 


o* 


Hi 


m 


-d 


CO 


o 




O 


o 


O 


fM 


.4 


Gt 






c 


nT 


W 


jj 


3 










• 










r>J 








3 


1 


i 










1 


•1 


♦ 


r<> 






o 


— 


O 




o 


c 


o 


O 


3 


o 




♦ 


fM 


in 


si- 


si 




o 




o 


o 


G 


G 


G 


.31 




O 


c 


03 






O 




1 








1 <J 






X 


r- 


O' 


n 


3 


G 




1 








1 


GI 


r\j 








• 


• 














in 


di 


fM 


















1 


1 






vn 




nT 


o 


>r 






o 






j 










3 




X 


3 


X 




M 


3 






1 


[ 








X 


















1 












X 












3 




























in 






! 










3 














vn 






1 






































3 
















1 
















03 


•4J 


1 r- 


fM 


IM 


w 




•si 


O' 


fU 


sT 


UJ 


o 






X 


Ul 


o 


o 


M 






U 1 


vTi 


si 






G 




CD 


o 


'O 


r- 


O' 


o 


ac 




c 


fSf 


O' 


fM . 


r» 


^ i 


1 




fM 


.r\ 


^r 


•d 


o ■ 


ar 


O' 


«; 


rd 


O' 


n 


o; 




« 


r'- 


<-* 


O 


lil 


o- 


<1 


fM 


c 


sT 


o 


O' 1 


o 


dy * 




>r 


O' 


<o 


cr 


o 


^ i 


< 


n 


cr' 


O 


IM 


n 


oi 








. ® 


O 


o 


o. 


3 


O 




f»S 


-H 




3 


k— 


—t 


3 


X 


■ n 




O' 


3 1 


X 


—4 


O 




si 


O' 


oi 






h- 


. »n 


— 


O' 


o 




>f 


— 


rfi 


U3 


O' 


o 


< 




3 


m 


X 


—4 


fd 


O 1 




O' 


d. 




n 


X 


ol 


X 
















fM 


/M 


Hx 


3 


fM 1 


m 


U 1 


X 
















•s. 


Ol 


*3 


3 


O' 


31 




o 


'T 


1 fM 




o 




'T 


fM 


rr 


3 


>r 


3 1 


6 


JJ 




U'3 


O 


fM 


—4 


o 


3 


si 




m 


fM 


fd 


n 


Ol 




fM 




1 o 






o( 
















ac 




— • 




dd 






O' 




• 


• 


• 


• 


• 
















'T' 




*M 




X 




fM 


o 


X 1 








1 






n 


X 


r- 


o| 


3 




3 


o 






i 


1 










-r 




O' 






o 




















3 


1 








3 




vJ" 






sT 




C-y 












vn 


<sJ 1 


03 


<-» 


— • 


CO 


O' 


fj 


O ' 






j 








vn 




■T 


f»1 


1 «-> 


fM 


•1 


3 
















JT i 


3 


O 


X 


rsj 


UN 


M 


3 , 






















1 — < 




13 


O 






1 










3 




sT 


X 


sj- 


o 


X 


o 






I 










* 




m 




<v 


o 


o 


















* 


fM 


X 


hd 


si 


JS 


3 1 






1 














o 


o 




o 


o 




















X 


n 


X 


_4 


*d 


3 ) 
















—* 


















o 


rri 


fTi 


cr . 


3 






• 


• 




4 


4 


• 1 




fM 


oi 


.4 


X 


X 


3 






'll 


rr» 




— i 


t3 






— 


r 




O ' 


3 






>r 


o 


rvt 


— 


3 


3 


a 


o 


•h 


vn 


fM 


o 


3| 








i 






o 


-0 


>r 


r- 


sf 




3 . 


3 




X 












3 ! 


< 


si 


fdl 


X 


3 


T) 


“>! 




—t 


1 








n 


3 




T! 


O 


3 




3 










j 






n 


3 


X 


3 


0 


3 


IM 


O' 


























>r ' 


O 




















si 


m 




X 


_4 






O 




m 


•J- 


— < 


-> 












> ■ 


3 






3 


X 


•X 


un 


■T 


3 




n 


si 


-i 


fM 




ol 




o 






O' 


fM 


o 


rf1 


U' 




<3 




O ' 


O 






fM 


o 


3 


si 


O' 


3 


3 


—4 


si 


o 


X 


T 


o, 


-r 


o 


oa 


r- 


/\ 


(M 


o 


















'T 


O 


sT 


m 


<n 




O 




• 


•1 








4 




• 


vn 




IM 


—4 


C ii 




— 


m 


t 






U> 






• 


in 


vT 


c 


o 


o ( 


X 


X 


O' 


o 


d4 


fM 


3 


v/> 


^r 




• 


■ 


• 


• 




“M 




3 






O 




03 


— 


« 


' • 


• 


• 










3 






o 






'll 


o 


• o 


o 








i 








in 






3 




si- 


o 


o 


t 














in 




<•1 




1 f'j 


























3 




3 






o 




















1 


i f^' 






n* 

1 






1 






1 








3 




1 






o 




















1 

o( 


fA 


'T 




1 

O 






o 


't 


■n 




O 






t 


sT 


1 

k 


O' 


in 


3 1 




m 


xl 


o 


O' 


n 


o 




C' 


f-i 


r 


O' 


M 


o» 






3 


•r 


si 


o - 


3 






X 


fM 




X 




^ ! 




fd 


ini 


-r 


X 


si 


3 


CM 


O) 


ai r<i 


'J' 


LT 


r» 


< 




r- 


n 






J 




fM 


3 




3 


3 


IM 




< 


(M 




.3 


X 


•n 


3 




• 


ca 




_4 


o 


3 




3 


tr 




m 


n O 






• 


O 


m 


fM 


o 


o 1 


X 


f'J 


s3l 


X 


n 


iM 


3< 


<y) 








• 




• 






in 


in 


o 




C3 




03 


m 


• 












X 


f'l 


in 


n 


•d 


Ol 




r\ 




O 


o 


o 






r 


c 


n 




X 1 


3 






o 


-d 


3 


3 


o 


1 




O 


Xl 


O 


rd 


•d 


3( 




o 




»M 






o 


M 


03 


-« 


si 




r~ 1 


O 


1 




o 










o 


Ti 


X 


3' 


vn 


-i 


G 


Ol 








fJ 






n 
















j 














n 






























j: 




m 




.-i* 


o 










1 








X 


3 


ot 


X 




M 


Ol 


















cr 


















1 




















3| 










•« 


03 


o 












in 






fM 


—• 


fM 


3 


X 


C [ 














m 






n 


o 




-< 


Oi 




















O' 


3 


X 


in 


—* 


3 
















^r 


—* 


u 






<-i 


















1 


r 




O 




O 


CJ 


o ! 






1 












• 




f- 


IM 




CJ 




















« 


'i* 


r- 


'J 


3 


O i 
















y< 






























X 


n 






• 


• 


4 


















o 




fM 


v_; 


O 


o 




''J 


ca 






m ra 






o 




fM 


o 


'X 


J , 




X 


xj 


O' 


in 




o 




fn 






1 




o 


j: 


•• 


r* 




-o 


*3 1 


3 


1 




■ f3 




3 


1 




3 i 


■ C 


3 




n 


'O 


— 


■G 




rO 


i 








in 


< 


i*j 


O' 








3 


1 




3 


1 


j 






IT j 


-t 


M 


rd. 


fVI 


■i 


si 


O 
















S3 


o 


in 


_ 4 




cr\ . 


3 








1 


1 








X 


o 




o* 


vn 




o 


















o 


'» 




in 


r 


3 


j 






1 












_d 


nj 


n 


fM 


si 


31 












<o 


(.j 




o 


rj 


r 


Al 


o 






•s. 


'T 


ru 




4> 


O i 




n 




n 


3 


3 






o 


IM 


o> 


O 


03 


3 




tD 


n 


in 


in 


in 1 


U 


j 




o 


fM 


U3 


X 


X 


O 


3 ^ 


3 


is/- 


M 


si 


O' 


g| 


rsl 


o 


r 


o 


IM 


fH, 


O 


















fM 


o 


si* 


»") 


'• J 


•h 


o 




• 


4. 


• 


• 


4 


41 




• 


o 




f 3 


_ 








fM 


-3 


—4 


.fM I 


3 






• 




3 


o 


3 


o 


t 


o 


3( 


X 


3 


3 




X 




« 


t 




• 


« 




■o 




m 






3 




X 


3 


• 


■ • 




• 
















O' 




es 






o 


o 


«-i 












in 






X 


—4 


3 


o 


f3 


G , 














in 




Hi 




rsj 


1 




r> 




















o 


nj 


1 




3 . 






1 






















n 




















'*• 










in , 




























Ui 












t 


.J 


1 






1 








'u; 














UJ 








« 


f/1 




•si 








• 


«/> 














• 


03 




•si 








« 


n 




“si 


a 






> 


1/1 






o. 






> 


03 


• mm 




oc 






> 


•/3 


• 


3 


■X 






> 


Ol 


• 


31 








OJ 


U.' 












.o 


sk- 


<n yi 








Z 1 


LU 


Oi 


O' 


O' 


>— 






Uj 


lU 


vn 


O'! 


r 




< 










r 






a 




3 






z 




< 1 


o 


Z 


C3 




r 






G 


Z 


n 


a: 




UJ 




C 






ac 




•li 




< 




-J 








Ui ‘ 










iC 








I 




-4*1 


•Q 


<1 


a. 


• 


Vi. 


a 




< 


< 


3. 


• 




or. 1 


a. 


1 


< 


< 


2. 1 


• 




2C 


a. 


< 


< 


r 


• 




i: 




o. 


Oj 


w 




X 




a. 


a. 


Sli 






X 


3 


2. 




i. 


UJ 


— 1 


H- 




D 


2. 


a 


X 




3 




3 


>i 




X. 


•n 




/I 




W3 




0. 


or 


1/1 


03 


03 






2. 


03 


03 


•03 




/I 








O' 


Ol 


X 


Oil 

1 



134 



jfiguie 26d QUEUc Wim EAK I A UTCREGEESSI VE SERVICE ‘lltES AND fCISSCN INtUI. TAEUIATICN OF SAtFlE STAIISTICS F08 ISt 
LISlBIrUIICNS OF CLtULAIED 1 NTE 5 AEfil V AL TiaES AT N = 5COO AND N= 10000 (I 2 ANC X 4), CUMULATED SERVICE TIMES (3 2 AND S 4), 
WAITING TIMES, WITH CR WITHOUT ZEROS (WI ♦ 0, WI, I FOE CASE N=2500, 5000, 1500, lOCOO) , CUMULATED WAITING TIMES (W1 BAR €tc.) 
IhL THE AVERAGED DIEEEREbCES BETWEEN S2 AND X2 (D2); OBTAINED FRCM ai = 500 BftilCATlCNS OF THE RUN SiSS#S8; RX==2.S7; RS=3 . 



2. As foi the steady state of the value cf E[ W ] and Z[W] 
a large N was required for high t and/or r, also for the 



steady state of the distibution cf W and a large value 



of N was required. Informally testing their distributions, 
ty using the Plotting Subroutines and the values of their 
parameters skewness and kurtosis, we may conclude that: 



(i) . The distri 
exponential form, if t 
under EXFLT and the 
about 2 and 6 respecti 
true only if the see 
net evident i.hat their 
has not been reached 

where the W has n 

ICOOO 

can obtain nothing abo 



bution of W, given that W>0, has an 
he plot results in a straight line 
'/^and 'jj 2 _ paramerers have the values of 
vely (see figures 23a-23d) . This is 
ady state has been reached while it is 
distribution is if the steady state 
. Thus observing the figures 24a-24d, 

ot reached the steady state yet, we 

ut rhe distribution of S. 



Note that 
distribution of 
distr ibution 
was one cf t 



for the H/M/ 1 
W, given that 
t known for the 
ectives of this 



is no 
he obi 



queue the steady 
W>0 , is exponential 
correlated queue and 
srud y . 



state 

The 

this 



(ii) . The distribution of T7 has the normal form if we 

get a straight line under NOBMPL and if the values of 

skewness and kurtosis are around 'O'. The nornality 

assumption holds for the steady state only, while the 

non-steady state distribution looks like an exponanxicnal. 

figures 27a-27c give us the feeling that has an 

. 10000 

exponential form, since a straight line appears under EXPLT 

(see figure 27b) and the parameters have the values 

2.3 and 7.3 respectively. These are close to the actual 
values 2 and 6 of the exponential distribution. We can see 
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from th€ flcts that as the T7 gets close to the steady 

state, it leaves rhe exponentional form and comes closer and 

closer tc the normal form. Figures 28a-28c, 23d shew us 

that U has already reached the steady state and 

10000 

furthermore its distrihution has eventually taken the normal 

form, since a straight line appears under NORMPL (figure 

28c) and the skewness and kurtosis have the value .7 and .8 

respectively. Compare its distribution with the 

distribution of 1? of figure 27b where U is not in the 

1 0000 

steady stats. Again we may state here that convergence to 

the steady state of the distribution of W and ¥ requires 

much larger N than convergence to the steady state cf their 
mean values for the same t and r. And also the higher t 
and/or r the larger N should be. As an example we can see 

that even though the steady state (in value) of T! for 

S#25#98 type has been reached fer N=10000, the steady state 
in distribution has net been reached yet. See figures 
24e-24g where we can observe that rhe disrribution of "R 
starts to leave the exponential form and to go tc the 
normal form. 
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Figure 2 
AND P0I5 
WAITING 



7a - QUEUE 
SON INPUT. 
TIMES ¥ 



WITH EAE1 
HISTOGRAM 
FROM 



20 0 0 



AUTOREGRESSIVE SERVICE TIMES 
OF THE CUMULATED AND AVERAGED 
THE RUN S#50#90; m=500 



REPLICATIONS, R X= . 2 ; RS=. 4 . 
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NUUbEK Cl OaUEHEO EAIHS > 500 



crj tn 
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Figure 27t - QUEUE 
AND EOISSON INPUT. 
CUMULATED AND 



fcITH EAR1 AUTOREGRESSIVE SERVICE TIMES 
EXPONENTIAL PLOT (EXPLT) OF THE 



AVERAGED WAITING TIMES 
S#5C#90; m=500 REPLICATIONS, RX=.2; RS=.4 



H FROM TEE RUN 

2 0 00 
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Figure 27c - QUEUE 
AND POISSON INPUT. 
AND AVEBAGZC WAITING 
HEPLICATICNS , BX=. 2; 



WITH EAE1 AUTOBEGBESSIV 

NOBMAL HOT (NOBMPL) OF 

TIMES TI FBOU THE BUN 

2000 

RS=. 4 . 



E SEBVICE TIMES 
THE CUMULATED 
S#50#90; m=500 
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Figure 28a - 
AND POISSON INPUT. 
WAITING TIMES H 



QUEUE WITH EAE1 AUTOREGRESSIVE SERVICE TIMES 
HISTOGRAM OF THE CUMULATED AND AVERAGED 



20 00 

REPLICATIONS, RX=.2; RS=.4 



FROM 



THE 



RUN 



S#50#25; m=500 
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NUtiBKR Of OROCBEO SAIRS » 500 



Figure 28t - 
AND POISSON 



QUEUE 

INPUT. 



«ITH EAR1 AUTOREGRESSIVE SERVICE TIMES 



EXPONENTIAL PICT (EXPLT) OF 



TH] 



CUMULATED AND AVERAGED WAITING TIMES 
S#5C»25; m=5CC REPLICATIONS, RX=.2; RS=.4 



^ FROM THE RUN 

a 0 00 
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Figure 28c - 


QUEUE 


WITH EAR1 AUTOREGRESSIVE SERVICE TIKES 


AtiD 


ECISSCN 


INPUT. 


NCRMAL HOT (NOEilPL) 


OF THE CUMULATED 


AND 


AVERAGE! 


WAITING 


TIMES U FROM THE 


RUN S?50#25; ni = 500 








2000 




RiPLICATICNS 


, RX=. 2; 


RS = .4 . 
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Centinuing 


the analysis of the results we present 


teloa the re 


isults we 


got 


from all runs. 


concerning the 


values of £[1!] for various 


values of t, r. 


N. 


N 


t 


r 


E[TI ] ( correlated) 


E[TJ](a/M/1) 


2000 


.25 


.25 


.3720 


.3332 


2000 


.25 


.50 


.4326 


.3332 


2000 


.25 


.90 


.7578 


.3332 


10000 


.25 


.95 


1.0576 


.3332 


10000 


.25 


.98 


1.7077 


.3332 


2000 


.50 


.25 


1.2106 


.9985 


2000 


.50 


.50 


1.581 1 


.9985 


2000 


.50 


.90 


5.0785 


.9985 


10000 


.50 


.95 


9.3196 


.9985 


10000 


.50 


.98 


21.0542=*' 


.9985 


10000 


.95 


.25 


23.5651 


18.3528 


10000 


.95 


.50 


33.3367 


16.3528 


10000 


.95 


.90 


123.9350* 


18.3528 


10000 


.95 


.95 


196.2460* 


18.3528 


10000 


.95 


.98 


342.9240* 


18.3528 


10000 


.99 


.25 


62.3285* 


55.3371 


10000 


.99 


.50 


77.9370* 


55.3371 


10000 


.99 


.90 


193.91 13* 


55.3371 


10000 


.99 


.95 


271. 7713* 


55.3371 


10000 


.99 


.98 


426.6750* 


53.3371 


Note that 


these V 


alues 


have been scale 


d from the values 


ettained in 


the si HU 


lations 


by multiplying 


by RS. Thus for 


t=0.50 and 


r=0.25 


the table gives E[ U ] (correlated) as 


1.2106; the 


value in 


Figure 


23d is 3.026x0 


.4=1.2106 . The 



values mariced * have net converged. 
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k study on these results gives us the following 
observations : 

1. For given traffio intensity t the value of E[U] is not 
oonstant, as it is for M/M/1 queues, but it depends on the 
value of the oorrelation parameter r. Thus we see that the 
value of E[U ] increase, as the value of r becomes larger 
and larger. Furthermore from figure 29a, where E[U] is 
plotted versus r we can see that the rate of increment is 
not constant (therefore is not a linear function of r) but 
it increases with r and eventually goes up infinitely as r 
goes to 1. 



2. For given correlation r, we see that the value of FCiJ] 
increases with t, rhe same as happens in the U/M/1 queue 

(see figure 29b). Eut regardless of the value of r, E[ « ] for 
M/a/1 queue is always less than for the EARMA model. Thus 
we may state that the weiiting time of the aurocorr elated 
service time model is greater than the corresponding waiting 
time of the M/M/1 model. 
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Figure 29a - QUEUE WITH EAR1 AUI0EZGEES3IVE SEEVICE 
IIHES ARC POISSON INPUT. PLOT OF E[ ¥ ] VERSOS r. 
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E[W] 



t=.98 




figure 29b - QUEJE WITH EAH1 AUTOfifGSZSSIVE SSFyiCS 
TiaES ANE POISSON INPUT, PLOT OF Ef TJ ] VERSUS t. 
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3. The values of £[¥] for the type S#99#98 model 

increases with N following the curve in figure 30 a fact 
that is known to happen in the a/a/1 queue also. We can see 
from that figure that the rate of increment comes dcwr as N 
increases and hopefully should come to zerc as ¥ reaches 

the steady state. Ccmparing then the M/M/l queue with the 

correlated queue we can see that the M/M/1 queue ccmierges 

faster than the correlated especially when r is high. Thus 

we can see frcm figures 32a-32c that both the value and 

distribution of the 7 for M/M/1 queue have reached the 

10000 

steady state, but the corresponding U cf the correlated 

10000 

queue (figures 33a-d) has not (at least in its 

distribution) . (See the skewness and kurtcsis values for W1 

EAR, W2 EAR, W3 BAR and W4 BAR in Figure 33d. These appear 

to be decreasing tc zero, but have certainly not reached 

zero by fl=10,000.) Looking at figure 31c, we can obtain 

the result that neither U for the correlated queue, or 

10000 

1! for the unccrrelated queue has reached equilibrium. 

10000 

This is because the skewness and kurtcsis values are 
definitely net close tc zero. 
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figure 30 
AND POISSON 



QUEUE WITH EAR1 AUTOREGRESSIVE SERVICE 
INPUT. PlOr OF E[17] FROa THE EON S#99#S8, 



300000 



TINES 
vs N . 
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Figure 31a - QUEUE 
AND POISSON INPUT. 
WAITING TIKES H 



KITH EAE1 AUTOREGRESSIVE SERVICE TIMES 
HISTOGRAM OF THE CUMULATED AND AVERAGED 
FROM THE RUN S*95#95; ai=500 



1 00 00 



REPLICATIONS, RX=3.8; RS=4 
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Figure 31b - 
AND POISSON 
CUMULATED AND 



QUEUE MIIH EAR1 AUTOREGRESSIVE SERVICE TIMES 
INPUT. EXPONENTIAL PLOT (EXPLT) OF THE 
AVERAGED WAITING TIMES ¥ FROM THE RUN 

l 0000 



S#95#95; m=500 REPLICATIONS, RX=3.8; RS=4 
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Figure 31c - QUEUE WITH EAR1 AUTOREGRESSIVE SERVICE TIMES 
AND POISSON INPUT. NORMAL PLOT (NORMPL) OF THE CUMULATED 



AND AVERAGED WAITING TIMES "W 



l 0 00 0 



FROM THE RUN S#95#95; 



m=500 REPLICATIONS, RX=3.8; RS=4 
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Figure 32a - a/M/1 QUEUE, 

AVERAGED y^AITING TIMES 'i 

REPLICATIONS, RX=4 ; RS=8 . 
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Figure 32b - M/M/1 QUEUE. EXPONENTIAL 
THE CUMULATED AND AVERAGED WAITING TIMES ¥ 
S#50#95; in = 500 REPLICATIONS, RX=4; RS=8 . 
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Figure 33a - QUEUE 
AND POISSON INPUT. 
WAITING TIKES ¥ 



1 0 0 00 

REPLICATIONS, RX=4 ; RS=8 



WITH EAR1 AUTOREGRESSIVE SERVICE TIMES 
HISTOGRAM OF THE CUMULATED AND AVERAGED 
FROM THE RUN S#50#95; in=500 
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Figure 33b - QUEUE 
AND POISSON INPUT. 
CUMULATED AMD 



WITH EAE1 AUTOREGRESSIVE SERVICE TIMES 
EXPONENTIAL. PLOT (EXPLT) OF THE 



AV ERAGED 

S#50#95; ni = 500 REPLICATIONS, RX = 4; RS = 8 



WAITING TIMES ¥ FROM THE RUN 
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Figure 33c - QUEUE WITH EAR1 AUTOREGRESSIVE SERVICE TIMES 



AND POISSON INPUT. 



NORMAL PLOT (NORMPL) OF THE CUMULATED 
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10000 



FROM THE RUN S#50#95; 



158 



-5CALE : = 0.960E-01 UNITS 



CO»H£LAThO CU£^U£ 



»“ j3 ^ O O 



(N O ^ 0 — 



o r* <7* 

<7' fS vO 

in <N r- 



•-* o cs 



® in 



o =7» o 



(n CN 



wi 

o 



O 'A ■» <N <0 

N in ^ ^ 

^ in ^ ;?> 

n »- in -o o 

<N o -N ■n o 

0 0 0-^0 



(M « •- m o 

»- f*' o <N 

B O r» -N 

^ rvi r* o 

O 35 <n 

O O O <N <n 



infsoNin*- 

in »- p- rvi 

CD in o 

aj rM in o 

ri| o rn »- 



o 

o 

o 

-5 

o 

.N 

o 

o 

o 

o 

o 

o 

o 

in 

o 

o 

o 

o 

o 



o o o 



m ~o 



in 



wi 

o 






CS sO 

in 

rn O' 

vD 

o -- 
-- o 



>o in 

UO <7^ 
if o 
rsi <N 

m o 



<7^ o n o 

BOBO 
B *- O O 

B B 3 O 

■B o r- o 

fS rn 'N O 

•- o 
in 

^ <N o 

B B O o 

rn in B cj 

cr» B »- o 

(N in »- o 

CN B B 

it 



yi B B in O 

O rn <N ^ O r-> B 

in in B r* »- o 

♦ ^ r- »- p' o 

^ in PM o 



»- o fN rn r« o 

at f- o 

m 



»- CN B r- cj 

rM »- rii B o o 

m D> in <71 o o 

j> — B in o 

^ r^l :» B PM O 



fM o rii *- 

n 

r^l 



B 

O 






tfi O' 

»- o 

o> ^ <y> 

O <N O 




^ O 



'-'NT' 
O B — 

B o* 

. o 

rN • • 

o PI -o 

• <N r- 



r» P5 o 

o p- o 

B B O 

O O O 

o -n 3 



;T 'N O 

<N O 

in 

TH -B O 

P- O' O 

P^ B o 

o o 

o o o 

t ca 
n 



B o 

O P- 

<-( JT 

B P5 

o o 
o o 

0 o 

d <d 

1 



B B 
B B 

in <N 
P- B 
f*' fN 
O O' 
o o 

o o 



o o 

B O 

B O 

P>l o 
O' o 
o 

0 o 

d d 

1 O 

in 



O' in 

?N *- B 

»- O O 

r>i in "N 

— o — 

o o »- 

o o o 

0 o o 

1 



JN 

B 

O' 

B 

B 



O 




O' o 
f-> o 
»- o 



0 o 

1 o 
m 



B B BOBO 
O' n O' O' O' o 

B in o o 

O' B o B B O 

B B P'- r* B o 

B — B B B O 

— o -n 0> n 3 

— o o o o 

o 

in 



3 -N O' n O' o 

B B n B P- o 

*- o n o 

B p^ -N n »- o 

B B B B O O 

B — 3 Ol O 

— O B O B 3 

— O O — 3 

O 

n 



bj 

3 

CU 

o 

o 



<N m p- r- 

B ® B O' 

n -- B O' 

»- B B 

O PM o 

0 0 0 -5 



n o 
p** o 
m o 

P5 O 

in o 

d o 
«- o 
n 



(N B B B B o 

Oi B P- B n o 

i^ fN B B <-? 

f'< o O' B o 

CN 3 N -n o 3 

O O O »- »- P'1 

n 

pij 



CO 

o 



p't 



3 



3 






n P- B PM B 

-5 »- n rn m 

«N oj PM -n n 

3 B O' 

— O ^ B 



o 



O O O 



I-* B vB 

“N O 

--BO 
<n »- rn 

-N O "N 

d o o 



B PI B 

cri — 

— > O' B 

PM O 3 

»- 3 -M 

O O O 



PI P- B 

P- r* B 

n p- — 

• n iP 

® 

B O PM 



ni B 3 

O 

n 

m p- w 

P- P5 o 

r- B o 

3 B o 

D B 3 



-M 



PVI P5 
<n B 
O' — 
P' O' 
PM B 

*M n 



B B 
B O' 
3 PM 
O O 

3 O 



B r» o n O' o 

B P^ n P5 o o 

n PM B B B O 

0 r* P- n P^ o 

3 3 PM O' CD 3 

O O O B 3 O 

O O O "N O O 

3 0 3 3 3 0 

1 I O 

n 



PM 

O 



»- O' B o 

r** -5 O' p» n 3 

■nr- B p- n o 

3 — B B B 3 

O O O P- P- o 

0 o o O' n o 

3 o o o o 

0 0 0 0 3 0 

1 I o 

n 



PM ^ O B O 

as B O P- PM B o 

•c m -n PM o 

a B -5 O' -5 O' o 

B 3 B P' P' 3 

■N O O O 3 

B - 00-500 

3 0 0 0 0 0 3 

o 

n 



— — B n CN o 

as B PM PM B 05 o 

-5 ^ r- B r- o 

3 B “5 B "N 3 O 

B 3 P- — O' O 

PM 3 O 3 O O 

P-» — 3 O — — O 

3 O O O O O O 

O 

n 



PM 



W1 



B O P- 

B r» B 

B — n 

. o o 

B t • 

— PM B 

o n 



O' o 
BOO 
3 W O 
0 — 3 

O O O 
I O 

in 



B B 05 O 

— n — o 

O' 3 3 

-M *M — 3 

B* d o o 

PM I 3 

n 



B O' O O 

PM -5 B B 

n p- o *M 

• p- -n o 

O' • • • 

B O P- O 

■N — I 



?N O 
CD 3 
r- 3 

— o 

3 3 

3 

. i/^ 



® — 

n n 

-5 r- 
B P- 
B PM 
PM PM 
— 3 



O 



3 — 

— B 
r- B 
•M — 
O' B 

o n 

n B 



3 



O 3 
3 O 
B 3 
B O 
-5 O 
-N 3 
B 3 



3 O 
3 

n 




!N 



W1 




n 

PM 

CD 



B 



3 -5 o 

P- P- 3 

PO (M 3 

PM O O 

0 3 0 

O 

n 




PM 



P- PM B — Pp O 

B -5 3 O 3 O 

3 — 3 r- BO 

PM B PM B B O 

B 3 O' fp* n o 

PM 3 O 3 3 

— 30330 

d d o d d o 

o 
n 



B — 

n 3 

O B 
PM n 
B O 
PM 3 
— O 

3 d 



n B o o 

tN n M o 

3 0 0 0 

3 3 3 0 

PM O 3 3 

- 3-0 
O B 'B O 

0 0 3 0 

O 

in 



33 

H 

C 

as 

<4 

Oi 




< Q 



as 5j 

< n • 

M H 

u cn cn 



cn 

<n • 

uM iO 

as o 

3 H 
lu as 
se 3 
cn as 



OU 

M 



cn 



03 

H 

TJ 






. 'B Nj 

^ » Wl . 

sc B uJ cn B 

-c a 2 J 3 

'B -■ H 

. SL4 -a u. 

-- H a: 3 n 

B B B as B 



B 

ca 



04 



OU 

. B P4 

^ ^ yj . B 

2 UJ B B B 

< Q 75 3 

3J at H B 

c • u iB U 4 

p- H as 3 i: 

B B B a: B 






B 



H 

B 



B 

3 



159 



figure 33<i - QUEUE kllh EAR1 A U'XCBtGRESSI VE SERVICE lltES AN! ECISSCK INEUI. TAfcULATICN OP SAKFIE STATISTICS POR THE 
LISIEIEUTICAS Ot CliCUlATED INTEtARRIVAL TI^ES AT M = 5000 AND lOOOO (X 2 AND X 4), CUMULATED SERVICE TIMES (S 2 AND S 4), 
WAITING TIMES, WITH OR WITHOUT ZEROS (UX ♦ 0, WI, 1 POE CASE N = 250C, 50C0, “iSOC, lOCCO) , CUttUiATED kAITIKC TIMES (Wl EAR etc.) 
AND THE AVERAGED DIEEEREACES BETWEEN S2 AND X2 (D2) ; OBTAINED PRCfl b = 500 EEtHCATICNS OF THE BUN S#50i95; RX=4; fiS=8 . 



Ancther result obtained from the simulation 

concerns the comparison of convergence of W and ¥. From the 

programs we have run we observe that 17 converges slightly 

faster than W (at least for low t and r) . Also the standard 

deviation of TJ is much less than of W. In figure 24d we can 

see that 'H has already reached the steady state but not 

10000 

li . Compare also the values of their standard 

10000 

deviations (.005 versus .2). Figure 34 gives the results of 
the type S»50#98 model for N=10000. We can see that while 

H is already close to the steady state, the fact that W has 

not reached it yet forces us to continue the simulation by 
increasing the value of N. 
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The model for the cross-correlated queue is: 



{X ,} = g , 

1 1 



jSJ is EASaA(1, 


1) over { 


E , EX/RS 


1 






i i 


That is; 




w . p . 




s. = \ 


1 ^ 






^ 1 


^ + 


A w . p 


0-6) 




1 


i 




where 


r rA 




w . p. r 




i-1 






A = ^ 








i 


^ rA 


+ 6 ax/Hs 


1 

• 

a, 

• 

> 




i- 1 


i 



RX/SS, 

i- 1 



} 



Again we chose the values .25, .50, .95, .99, for t 
and .25, .50, .90, .95, .98 for r and we run a toral cf 20 
runs to cover all ccibinaticns of t and r, naming each run 
as SC#tr#rr. 

Analyzing the results obtained from the simulation 
we observe the fcllowirg: 

1. The value of the traffic intensity t anc the 

correlation r is the main factor for convergence of W and U. 

That is, high traffic intensity and/or high correlation 

require N to be large in order for W and W to reach the 

steady stats. Thus simulating the models with t<.50 and 

r<.95 for S=10000 it was possible for H and U to reach the 

steady stats. In figures 35a-35d, 36a-36c we can see the 

justification of that conclusion. The values .1078 and 

.0957 of A , and 7 represent the steady state of 

10000 10000 



the SC»50#25 model. On the ether hand, regardless of the 
value of t the models with r>.98 did nor reach the steady 
state for N=10000 and regardless of the value of r the 
models with t>.95 did net reach the steady state for 
N=1C,000 . figure 37 gives us rhe non-convergence of the 
model SC#99#25. Because of the high traffic intensity and 
in order to get the steady state of that model we have to 
continue simulating by increasing the value of N. 
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Figure 35a - QUEUE WITH CHOSS-CORRELATEE SERVICE TIME 

SEQUENCE AND POISSON INPUT. HISTOGRAM OF THE WAITING TIMES 

W WITHOUT ZERCS FROM THE RUN SC#50#25; m=500 

X 0000 

REPLICATIONS, RX=2; R5=4 . 



No OF ZER0S=279. 
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figure 35d - QUEUE WITH CftOSS-COBUELAI ED SEBVICE TlflE SEQUENCE AbC POISSCN INPUT. TABUIATICN Of SAEFIE STATISTICS fOH T«£ 
£1 S15IEUTIC bS CJ? CUEULATED INTEI AFRI V AL TIdES AT N=5C00 AND N= 10000 (X 2 AND X 4) ^ CUMULATED SERVICE TIMES (S 2 AND S 4) ^ 
KAITING TIMES, UITh OR WITHOUT ZEROS (UI ♦ 0, HI, I FOB CASE N=2S0C, iOOO, iSOO, lOOOO), CUMULATED WAITING TIMES (HI BAB etc.j 
AND THE AVERAGED DIFFEBEbCES BETWEEN S2 AND X2 (D2) ; OBTAINED FBCM sn^SOO BEEllCATIONS OF THE BUN SCI50|25; EX=2 ; BS=4 . 



2. High trafic intensity and correlaticn affects the 
convergence of the distrituticn of W and U also. 
Furthermore the steady state of their distributions, 
reguires larger N than is required for the steady state of 
their expected values. In examining that ccnclusicn, look 
for example at the results of the simulated model SC#50#98, 
where although W and "J! have already reached the steady 

state for N= 10000, their distribution has not converged 

yet. The above results have been obtained from Figures 

30a-38c, 5Sa-39d where W and U are presented. We 

10000 10000 

can see from these figures (38a-38c) that although the value 

cf W has reached the steady state (see Figure 39d) , 

10000 

its distribution cannot be characterized since the plot 



under EXfLT (figure 38b) is not a straight line 
parameters away 
desired values 2 and 6 . The same results can be 
from the figures 39a-39c which deal with TJ 

10C00 



and the 
frcm the 
obtained 
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Figure 36a - QUEUE WITH 
SEQUENCE AND POISSON INPUT. 
AVERAGED WAITING TIMES ¥ 



10000 



CHOSS-CORRELATEE SERVICE TIME 
HISTOGRAM OF THE CUMULATED AND 
FROM THE RUN SC#50#25; in=500 



HEPLICATICNS , RX=2 ; RS=4 
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Figure 36t - QUEUE WITH 
SEQUENCE AND POISSON INPUT. 

THE CUMULATED AND AVERAGED WAITING TIMES ^ 
SC#50#25; m=500 REPLICATIONS, R X=2 ; RS=4 
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EillHATtD PARAfltTEBS OP DATA ; GArtflAl = 1.7bJ51SE-01 GAflBA2 = - 1 . 895754 E-Ol 



figure 36c - 



QUEOE 

SEQUENCE AND POISSON INPUT. 
CUMULATED AND 



WITH cross-correlate: service time 



NORMAL PLOT (NCRMPL) OF THE 



AVERAGED WAITING TIMES TJ FROM TEE RUN 

1 0 000 



SC#50#25; eb = 500 REPLICATIONS, RX=2; RS=4 
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LXSl£l£UlICl »5 CP CU£U 11 I£D I N I E £i Aafil V AI TIW£b AI N = SC 00 AND N= 10000 (X 2 AND X 4 ), CUMULATED SERVICE TIMES (3 2 AND 
aAlTlbG TIMES, UITH CR WITHOUT ZEROS (WI ♦ 0 , WI, I FOB CASE N^ 250 C, 5000 , T 500 , lOCOO) , CUMULATED UAITII^G TIMES (Ml DAB 
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Figure 38a - QUEUE WITH CROSS-CORRELATED SERVICE TIME 
SEQUENCE AND POISSON INPUT. HISTOGRAM OF IHE WAITING TIMES 



W WITHOUT ZEROS FROM THE RUN 

10000 

REPLICATIONS, RX=2; RS=4 . No OF ZEROS=245. 



SC#50»98; m=500 
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Figure 38b - QUEUE Wir H 
SEQUENCE AND POISSON INPUT. 



CEOSS-COHHELATED SERVICE TIME 
EXPONENTIAL PLOT (EXPLT) CF THE 

HUN SC#50#98; 



WAITING TIMES W WITHOUT ZEROS FROM THE 

100 00 

ni = 500 REPLICATIONS, RX=2; RS = 4 . No OF ZERCS=245. 
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fcSlIflATED fc>AbAntTEaS OF DATA ; GAMHAl = J.3445eSE 00 GAttflA2 » 1.099609b 01 



NOMUlin Ct ORDERED JAIKS = 255 



Figure 38c - QOEOE WITH CROSS-CORRELATED SERVICE TliSE 

SEQUENCE AND POISSON INPUT. NORMAL PLOT (NORMPL) CF THE 

WAITING TIMES W WITHOUT ZEROS FROM THE RUN SC#50#98; 

1 00 00 

ai = 500 REPLICATIONS, RX=2; RS = 4 . No OF ZERCS=245. 
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Figure 39a - QUEUE WITH 
SEQUENCE ANC POISSON INPUT. 
AVERAGED WAITING TIMES ^ 



CROSS-CORRELATED SERVICE TIME 
HISTOGRAM OF THE CUMULATED AND 
FROM THE SC#50#98; m=500 



1 0 000 



BEPLICATICNS , RX=2 ; RS=4 
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Figure 39t - QUEUE WITH CROSS-CORRELATED SERVICE TIME 
SEQUENCE AND POISSON INPUT. EXPONENTIAL PLOT (EXPIT) OF 
THE CUMULATED AND AVERAGED WAITING TIMES TJ FROM THE 



1000 0 



SC#50#98; a=500 REPLICATIONS, RX=2; RS=4 
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ES'IIHATED PARAMETEHS OF DATA : GAKflAl » 1.b0220tlE 00 GABKA2 = 14.912*109E 00 



NUnUtR eg OBDtftED IKIU3 « 500 




figure 39c - 
SEQUENCE AND 
CUMULATED AND 



QUEUE WITH CROSS-CORR 
POISSON INPUT. NORMA 
AVERAGED WAITING TIM 



ELATED 


SERVICE 


TIME 


L PLOT 


(NORMPL) 


CF THE 


SS w 


FROM 


THE 



10000 



SC#50#98; m=500 REPLICATIONS, RX=2; RS=4 . 
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figure 39d - QUEUE HITU CEOSS-COfiRELAT ED SERVICE TlttE SEQUENCE AKD POISSCN INPUT. TABUIAIICN Of SAEtlE STATISTICS fOD lUE 
EISTRIEUTIONS OF CUflUlATED INIERAHRIVAl HUES AT N==5COO AND K= 10000 (X 2 ANC X ^4 ) , CUMULATED SERVICE TIMES (S 2 AND S 4), 
WAITING TIMES, WITH Cfi WITHOUT ZEROS (WI ♦ 0, WI, I FOR CASE N=250C, 5000, T50C, lOCOO) , CUMULATED WAITING TIMES (W1 BAR etc.) 
AND THE AVERAGED DlffERENCES BETWEEN S2 AND X2 (D2) ; OBTAINED FRCtt 0=500 EEEIICAIICNS Cf THE RUN SC#50lSE; £X=1.S; KS=6 . 



3. To tast the distribution of H and ¥ we also usee the 



plotting subroutines and the values of their skewness and 
kurtosis parameters as we did with the gueue with dependent 
service . 



The results obtained from the 

fellows: 



analysis are as 



(i) . G 
exponential 
result in a 
around their 
the SC#50#25 
35a-35d) , wh 
straight 1 
fer and ^ 2 . 
expenen tial 
steady state 
still in a t 



iven that H>0, we may say 
form, since the plots unde 
straight line and the valu 



actual 


values 


2 and 


6 respect 


model 


( which 


is r 


epresented 


ere the 


steady 


state 


has been 


ine 


and 
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values 



respectively. Ke have to s 
distribution comes only if W 
, and we cannot have the same 
ransient state. 



that W h 
r EXPLT Subr 
es of /j,and 
ively. For e 
by the f 
reached, gav 
1.8 and 
tate here th 
has reache 
result if th 



as an 
cutine 
^^are 
xample 
igures 
€ us a 
4.7 
at the 
c the 
€ S is 



(ii) For the distribution of U we can say that the \ 
normality assumption is an appropriate one, since the 
graphs obtained under NORMEL result in a straight lire and 
the value of skewness and kurtosis is around 'O'. The 
results, for example, obtained from the SC#25#25 model gave 
us a straight line under NOEMFL and the values .17 and -.05 
fer the skewness and kurtosis respectively (see figures 
40a-40c) . The normality assumption holds if the ^ has 
reached the steady state. At the begining of the transient 
state it is not possible to characterize its distribution, 

while when U is in a transient state we can see that its 

2 

distribution takes an exponential or a X form. This is 
left little by little as U comes close to the steady state 
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and then eventually it takes the normal fcxm. The model 
SC#25#95 is a representative example of that result. In 



that model, although TJ has not reached the steady state yet 



(i!l=10000) we can observe the values of skewness having as 
3.5, 1.8, 1.6 1.5 and the values of kurtosis as 24.7, 7.2, 
5.7, 4.3 for N=2500, 5000, 750C, 10000 r espectivelly . That 
is the values start from a high level, then come down, they 
pass from the actual values of exponential 
model and then continue decreasing and hopefully in the 
steady state they should reach a value around '0‘ where the 
value of skewness and kurtosis for the normal distribution. 
Ihe same assumptions are obtained analyzing the plots of 
that model also. That is for N=2500 the plots are not 
linear under both EXPLT and NOEilPL, for N = 5000 a straight 
line appears under EXPIT, for N=7500 and 10000 the plot 
leaves gradually the straight line under EXPLI and begins to 
give us a straight line under NOBMPL. 
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Figure 40 

SEQUENCE 

AVERAGED 



a - QUEUE WIDH 
AND POISSON INPUT. 
WAITING TIMES 



CROSS-CORRELATED SERVICE TIME 
HISTOGRAM OF THE CUMULATED AND 
FROM THE SC#25#25; m=500 

0 000 



REPLICATIONS, RX=1.5; ES=6 . 
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Figure 40t - 
SEQUENCE ANC 
THE CUMULATEE 



QUEUE WITH CROSS-CORRELATED 
POISSON INPUT. EXPONENTIAL 
AND AVERAGED WAITING TIMES 13 



SERVICE TIME 
PLOT (EXPIT) OF 
FROM THE 

1 0000 



SC#25#25 ; 



tn = 500 REPLICATIONS , RX = 1.5; 



RS=6 , 



183 



' MMA1 OATA— r CAKMAl = 1 . 0 0 5 1 £ t 01 GAMMA? ^ . 9B£A I 0E-O2" 




l/> Lf. 




o o 
« 1 

^ uu 
<M q> 

ro ^ 
o 9 



It 



♦ «' 




t ( 
X 



Figure 40c - 
SEQUENCE AND 
CUMULATED AND 



QUEUE WITH CROSS-CORRELATED SERVICE 

POISSON INPUT. NORMAL PLOT (NORMPL) 

AVERAGED WAITING TIMES f! FROM 

1 000 0 



TIME 
CF THE 
THE 



SC#25#25; m=500 REPLICATIONS, RX=1.5; 



RS=6 . 
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Figure 4Cd - QUEUE Vllti CROSS-COFRELATEE SERVICE TIflE EEQUENCS AtD fOISSCN INfUl. TAbUlATICK OF SftFlE STATISTICS FOB 
tlSlf lEUTICNS OF CUrULAIED INTSBARRIVAL TIKES AT N = 5000 AND N= 10000 (X 2 AND X 4), CUHULATED SERVICE TIMES (S 2 AND 
WAITING TIMES, WITH Cfi WITHOUT ZERCS (WI + 0, WI, I FOR CASE M = 250C, 5000, 'iSOC, 10CC0), CUMULATED WAITING TIMES (W1 oAB 
AND THE AVERAGED DIFFERENCES bETWEEll S2 AND X2 (E2) : DETAINED FKCK b = 500 REfllCAIICRS CF THE RUN SCI25»25; RX = 1.5; ES=6 . 



4 . The 


following 


values for 1? have been obtained from the 


runs with 


various values of t 


and r and 


for 11 = 10000 


N 

1000C 


t 

.25 


r 

.25 


E[TI] 

. 1 169 


ErK](a/M/i) 

.3332 


10000 


.25 


.50 


. 1724 


.3332 


10000 


.25 


.90 


.3538 


.3332 


10000 


.25 


.95 


.4075 


.3332 


10000 


.25 


.98 


.4679 


.3332 


10000 


.50 


.25 


.3831 


.9985 


10000 


.50 


.50 


. 5648 


.9985 


10000 


.50 


.90 


1.4589 


.9985 


10000 


.50 


.95 


2. 1109 


.9985 


1000C 


.50 


.98 


3.7196* 


.9985 


10000 


.95 


.25 


10.4672 


16.3528 


10000 


.95 


.50 


13. 4784* 


16.3528 


10000 


.95 


.90 


29. 1312* 


16.3528 


10000 


.95 


.95 


75. 1180* 


16.3523 


10000 


.95 


.98 


145.3484* 


18.3528 


10000 


.99 


.25 


35.5986* 


53.3371 


10000 


.99 


.50 


43. 1616* 


53.3371 


10000 


.99 


.90 


98.9556* 


53.3371 


10000 


.99 


.95 


139.0167* 


53.3371 


10000 


.99 


.98 


219.7920* 


53.3371 


Analyzing the above lesults we may state conclusions 


as follows 


• 








(i) . 


For a given traffic 


intensity 


t the value of 1! is 


not constant, but depends on 


the value 


of correlation r, as 


happens in 


the queue 


with dependent service. We observe 


also that 


1 increases as 


r increases, but the rate of 


increment 


is larger 


as r goes 


to 1, and 


very small for low 
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values cf r; therefore TJ is not a linear function cf r. 
Figure 41 gives the flats of "R versus r for t = . 25 , .50, 
.S5, .98 and also for the a/M/1 queue (r=0) . Furthermcre we 
can observe that the value of 1? for the correlated queue, 
comparing it with the result for the M/M/1 queue, is not 
always greater as in the case cf dependent service, but 
depends on the value of r. It is less for r<.50 and greater 
for r>.90 

(ii) . For a given correlation the value of f increases 
as t increases. Furthermore high values of t cause high rate 
of increment, which becomes larger as the correlation 
increases. Figure 42 gives the plots cf T? versus t for 
r = 0, .25, .50, .90, .95, and .98 . Observing that 
figure we can see also that for a given t the order of 
increment of ¥ is for r = .25, .50, 0, .90, .95, .98 
respectively. Thus we can conclude that xhe waiting time 
fcr the cross-corr elat ed queue is generally greater than the 
waiting time for M/M/l queue when the correlation is high. 

5. The last result obtained from the simulation of that 
model was the convergence of H and TJ. (These statistics 
alsc converge in value and in distribution slower than the 
same statistics in a/M/1 queue). The mcdel SC#25#98 reveals 
that U converges faster than W (see figure 43) . But the 
model SC»S9#50 gave us the conclusion that W converges 
faster than "3 (see figure 44) . Thus we can state here that 
fcr low values of t converges faster than H and fcr high 
values of t the opposite is true. 
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Figure 41 - QUSUE WITH CROSS-CORRELATED SERVICE TIMES AND 
ECISSON INPUT. PLOT OF E[ IT ] VERSOS r. 
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Figure U2_ - QUEUE WITH CROSS-CORRELATED SERVICE TIMES AND 
POISSON INPUT. PLOT OF £[¥] VERSUS t. 
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figure 44 - QCEUE WITH CROSS-CORRELATED SERVICE TIrtE SEQUENCE AND ECISSCN INtUT. TABUIATICN OF SACELE STATISTICS FOR THE 
CIST BI EUIIONS Of CUMULATED INTERARRIVAL TIHES AT N-5000 AND N=10000 (X 2 AND X 4), CUMULATED SERVICE TIMES (S 2 AND 54)^ 
WAITING TIMES, WITH CB WITHOUT ZERCS (WI ♦ 0, WI, I FOR CASE N=250C, 5000, TSOC, 10C0C) , CUMULATED WAITILG TIKES (W1 BAB etc.) 
ANC THE AVERAGED DIFfEREfCES BETWEEN S2 AND X2 (D2) ; OBTAINED FRCK 0=500 EEEIICATICNS OF THE BON SC399I50; EX=2.97; RS=3 . 
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